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Abstract 



Turbulence is a major factor limiting the achievement of better toka- 
mak performance as it enhances the transport of particles, momentum 
and heat which hinders the foremost objective of tokamaks. Hence, 
understanding and possibly being able to control turbulence in toka- 
maks is of paramount importance, not to mention our intellectual cu- 
riosity of it. We take the first step by making measurements of tur- 
bulence using the 2D (8 radial x 4 poloidal channels) beam emis- 
sion spectroscopy (BES) system on the Mega Amp Spherical Toka- 
mak (MAST). Measured raw data are statistically processed, generat- 
ing spatio-temporal correlation functions to obtain the physical char- 
acteristics of the turbulence such as spatial and temporal correlation 
lengths as well as its motion. The reliability of statistical techniques 
employed in this work is examined by generating and utilizing syn- 
thetic 2D BES data. The apparent poloidal velocity of fluctuating 
density patterns is estimated using the cross-correlation time delay 
method. The experimental results indicate that the poloidal motion 
of fluctuating density patterns in the lab frame arises because the pat- 
terns are advected by the strong toroidal plasma flows while the pat- 
terns are aligned with the background magnetic fields which are not 
parallel to the flows. Furthermore, various time scales associated with 
the turbulence are calculated using statistically estimated spatial cor- 
relation lengths and correlation times of turbulence. We find that tur- 
bulence correlation time, the drift time associated with ion tempera- 
ture or density gradients, the ion streaming time along the magnetic 
field line and the magnetic drift time are comparable and possibly 
scale together suggesting that the turbulence, determined by the local 
equilibrium, is critically balanced. Finally, we argue that we have pro- 
duced a critical manifold in the experimentally obtained local equilib- 
rium parameter space separating dominant turbulent transport from 
a non-turbulent or weakly turbulent state. It shows that the inverse 
ion- temperature-gradient scale length is correlated inversely with q/e 
(safety factor/ inverse aspect ratio) and positively with the plasma ro- 
tational shear. Practically, this means that we can attain the stiffer ion- 
temperature-gradient, thus hotter plasma core, without increasing the 
rotational shear. 
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Chapter 1 



Fusion and turbulence in tokamaks 



We do not inherit the planet from our parents, we borrow it from our children. 

- Native American 

Imagining our world without any electrical power just for a few seconds, we 
come to a non-negotiable conclusion: we require electrical power to sustain and 
flourish in our lives. This power has been generated from natural resources, such 
as coal, oil, natural gas, and they inevitably produce carbon dioxide which may 
be argued as a cause of global warmingj^ Moreover, the limited amount of such 
resources is another vital concern [Ij. Thus, it will be undoubtedly beneficial 
if we can generate electrical power without further enhancing global warming 
and without depleting the limited amount of resources so that we can return the 
planet to our children without corrupting it further. Fusion power satisfies these 
criteria. 

1.1 Fusion in tokamaks 

Energy can be generated by fusing two nuclei of deuterium and tritium, 
which produces a harmless helium nucleus and a neutron with a total energy 
release of IT.QMeV ||2j. To fuse the two nuclei they must overcome the repulsive 
electrostatic force which can be achieved by heating them to (large thermal en- 
ergy) around ~ 10 — 20 keV at which reasonable performance of a fusion reaction 
can be achieved ||3l. Note that temperatures are expressed in the energy unit of 
eV throughout this work. At this high temperature, particles are ionized, i.e., 
become plasmaj^ 

ilPCC, Synthesis Report in IPCC AR4 SYR 2007 

^Being ionized is not a sufficient condition to be in the state of plasma. However, ionized 
gases in this work satisfy plasma criteria |^: i) Debye length is smaller than the system size; ii) 
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1.1. Fusion in tokamaks 




Figure 1.1: (a) Helical motion of a single charged particle along a magnetic field 
line in the absence of an electric field; (b) Cross-field motion of a single charged 
particle at the marked position due to cumulative effect of Coulomb collisions 
with other particles. Figure is taken from Ref. [|51 

1.1.1 Tokamak concept 

Given that such hot plasmas with a sufficiently large density are created, the 
performance of a fusion power plant depends on how long they can be confined 
within a finite spatial domain. One way to confine the plasma is using the Lorentz 
force, a basic concept of "magnetic confinement": 



where m, v and Ze are the mass, the velocity and the charge of the particle, re- 
spectively, and c is the speed of light. E and B are the electric and magnetic 
fields, respectively. By this law, a single charged particle in a strong magnetic 
field in the absence of an electric field is constrained to move along the magnetic 
field line with a helical trajectory, i.e., plasmas are confined in the perpendicular 
plane with respect to the magnetic field, as shown in Figure [LTj^a). 

If the confined plasmas were coUisionless in infinitely long parallel straight 
magnetic field lines, then a perfect confinement would be achieved. But, such 
field lines cannot be generated in practice, and furthermore, plasmas should not 
be coUisionless if the goal is to extract power from fusion reactions, i.e., particles 
need to collide with each other occasionally. 

Even though infinitely long straight magnetic field lines are not feasible, effec- 
tively "infinitely long" lines are certainly possible: a closed field line is infinitely 
long in the sense that starting and ending points are indistinguishable. This is the 
basis for a concept of the "TOKAMAK" a transliteration of a Russian acronym a 
toroidal chamber with axial magnetic fields. The geometrical definition of a torus, 
on which "infinitely long" magnetic field lines lie, i^ 

number of particles in a Debye sphere is large; and iii) typical plasma frequency is larger than 
collision frequency with neutrals. 
^From Wikipedia. 




(1.1) 
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Figure 1.2: Geometrical structure 
of a torus showing (a) poloidal 
(red arrow; short-way around) and 
toroidal (blue arrow; long-way 
around) directions; (b) a torus is a 
product of two circles: a red circle in 
a poloidal plane and a pink circle in 
a toroidal plane. Figures are taken 
from Wikipedia. 




Figure 1.3: Basic configuration of 
magnetic field lines in a tokamak: 
toroidal cage with closed field lines. 
Figure is taken from Ref . ||5j] 



a surface of revolution generated by a circle in three dimensional 
space about an axis coplanar with the circle 



as shown in Figure 1.2 A tokamak basically creates nested tori, where each torus 
is referred to as a flux surface, such that no magnetic field lines are connected be- 
tween the two tori within the Last Closed Flux Surface (LCFS) unless there exist 
radial magnetic perturbations. This means that magnetic fields can be described 
in two dimensional space, i.e., in poloidal and toroidaQ directions depicted with 



red and blue arrows, respectively, in Figure 1.2 'a). Although the poloidal cross- 
section of a torus is circular in its mathematical definition, that of a flux surface in 
a tokamak does not have to be a circle; in fact, it is usually 'D' shaped in practice. 
In addition, the flux surfaces are not usually concentric due to the Shafranov shift 
||6|, which means that the centres of the nested flux surfaces do not coincide. A 
circular-shaped flux surface in a tokamak, as an example, is illustrated in Figure 
O 



A coordinate system of a tokamak is illustrated in Figure 1.4 The axis about 
which poloidal cross-sections are revolved is defined as Z-axis, and its value is 
the height from the midplane (a plane containing the magnetic axes). Note that 
the centres of different flux surfaces are, to a good degree of precision, copla- 
nar and perpendicular to the Z-axis, i.e., the Shafranov shift is mostly in the in- 
direction where R denotes the major radius whose value is the distance from the 



^Toroidal fields are generated by external coils, known as the toroidal field coils, whereas 
poloidal fields are generated by plasma currents flowing in the toroidal direction. 
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Magnetic Axis 




Figure 1.4: A coordinate system of a tokamak showing the major radius R, minor 
radius r, height Z, toroidal (p and poloidal 9 angles, is a measure of plasma 
size at a flux surface ^. Figure is taken from Ref. [|5l 

Z-axis. The minor radius r is the distance from the magnetic axis, and is the 
half diameter of a flux surface at the magnetic axis height, i.e., at the midplane. 
Here, ijj is the flux surface label. We use a (without the subscript ip) to denote the 

of the LCFS, and Rq for R at the point where — )■ 0, which are measures of 
total plasma size. Note that r and can be different unless there is no Shafra- 
nov shift with circular poloidal cross-section. and 6 denote the toroidal and 
poloidal angles, respectively. The origin of (f) is not defined as it is not necessary 
for it being a symmetric direction; while 9 is measured from the outboard mid- 
plane. Outboard (inboard) is the region where R ■ r > {R ■ r < 0). Note that 
in some chapters of this work we use a Cartesian coordinate system with a local 
approximation, and it is explicitly stated when we do so. 

Collisions in a plasma are not like ordinary instantaneous collisions in the 
sense that they are due to the long range Coulomb force with stochastic inter- 
action between charged particles (binary) known as Coulomb collisions HZl- The 
characteristic Coulomb collision time of a plasma is usually considered as the 
time a particle takes to have a change of angle of order unity in the direction of 
velocity. Figure [LT| b) depicts the effect of Coulomb collisions where the marked 

5 



1.1. Fusion in tokamaks 



position corresponds to the cumulative effect of many Coulomb collisions. Al- 
though collisions are necessary for fusion power generation, it is obvious that 
they will degrade the confinement time. However, such a reduced confinement 
time due to Coulomb collisions is not what makes generating economical fusion 
power arduous since the associated particle diffusion is estimated to be at least an 
order of magnitude smaller than that experimentally observed Q- The observed 
transport is, in fact, anomalous exceeding both classical and neoclassical trans- 
port BSl by more than an order of magnitude, and it is believed to be associated 
with plasmas being turbulent in a tokamak 0, which has been the motivation of 
this work on plasma turbulence. 

As the final practical goal of the fusion community is to light up the whole 
world with economical fusion power plants, we must be able to ignite plasmas, 
i.e., generate self-sustained burning plasmas. Known as the Lawson criterion 
IITOl , the ignition condition can be written as a triple product of density n, tem- 
perature T and energy confinement time te of plasmas: 

hTte > 3 X 10^^ m-^ keV s, (1.2) 

meaning that have as high a density as possible, be as hot as possible and retain 
particles as long as possible. The International Thermonuclear Experimental Re- 
actor (ITEE[^ endeavors to satisfy the condition with n = 10^° m~^, T = 10 keV 
and Te > 3 s. 



1.1.2 Mega Amp Spherical Tokamak (MAST) 

The spherical tokamak is a type of tokamak whose main difference from a 
"conventional" tokamak is the aspect ratio e^^ = Rq/ a. A spherical tokamak has 
a tighter aspect ratio, Sq^ ^ 1.5 than a conventional one Eq^ 3 as illustrated in 
Figure 1.5 Once it was predicted that a small aspect ratio tokamak could achieve 
higher beta /3 operation [[TT1[T2| with a subsequent proposal of a new arrangement 
of magnetic coils (H in 1984, the Small Tight Aspect Ratio Tokamak (START) 
with ^ = 1.3 was built at Culham Center for Fusion Energy (CCFE), known as 
Culham Laboratory then, and the first plasma was reported in 1992 [14]. Here, /3 
is the ratio of the plasma pressure to the magnetic field energy density which can 
be interpreted as the higher the (3, the more economical it is because generating 
magnetic fields and building the coils are costly. Successful results from START 

^Step size of random walk is the Larmor radius due to Coulomb collisions. 
^Step size of random walk is the width of banana orbit of trapped particles due to the toroidal 
geometry. 

''www.iter.org 
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"^'Batge aspect ratio 
(eonvervtionaltokamak] 



Figure 1.5: Comparisons between spherical and conventional tokamaks. Geomet- 
rical configurations between the two tokamaks differ in their aspect ratio. Figure 
from CCFE. 



came out in 1997 that volume averaged (and central) /3 reached ~ 11.5 % (and 
~ 50.0 %) and that j3 could be increased further using higher input power [15] . 
The highest volume averaged (3 achieved, then, was 12.6 % held by DIII-E|^in San 
Diego, US ||16|. Subsequently, START had many experiments with the averaged 
/3 ~ 30 % and the highest (perhaps the highest record even to the current date) 
value of 40% [17]. 

These successful experiments on START led the team to upgrade START to 
the Mega Amp Spherical Tokamak (MAST) [[TSl with i?o ~ 0.8 m and a ~ 0.6 m 
(^Q ^ ~ 1.3) with plasma currents up to ~ 2.0 MA. The first plasma in MAST was 
reported in 2001 UM- Then, H-mode ||20l operation on MAST was reported ||2T|g| 
H-mode compared to L-mode (where H- and L- stand for high and low confine- 
ments) has a longer energy confinement time with steeper edge density and tem- 
perature gradients resulting in a larger stored energy in plasmas. In other words, 
H-mode plasmas establish an edge confinement "barrier" that reduces transport. 
H-mode operation has been achieved in many tokamaks (perhaps most of ex- 
isting tokamaks), and they are accompanied by a reduction of turbulence ([23]; 
and references therein). This, again, provides the motivation of this work on 
plasma turbulence: what are the mechanisms that suppress the turbulence and how do 
we achieve it^\ 

Before we sketch a simple picture of plasma turbulence, let us provide a pic- 
ture of MAST (Figure [L6) > and a table summarizing some of its parameters (Table 



1.1 1 from which the experimental data in this work were obtained. 



^DIII-D has i?o ~ l-7m and the central magnetic field 2.0T; whereas START had i?o = 0.3to 
with the central magnetic field ^ 0.5 T. 

'H-mode was also achieved on START II22I . 

^"H-mode operation is an empirical achievement without much physical understanding of it. 
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Figure 1.6: A picture of MAST. The vacuum vessel and other structures are CAD 
drawings, and the overlaid plasma inside is a real picture from a MAST dis- 
charge. Figure from CCFE. 



Table 1.1: Typical MAST parameters lEl 



Parameter 


Value 


Magnetic field on axis 


0.5 T 


Core temperature 


1.5 keV 


Core density 


5 X 10-*^^ m^^ 


Plasma current 


1.3 MA 


Pulse length 


0.5 s 


Aspect ratio 


1.3 



1.2 Plasma turbulence in tokamaks 

Our theoretical level of understanding of turbulence, or turbulent flow, is so 
minimal (but interesting^ that even the definition of turbulence is not well estab- 
lished ||25|. Consequently, experimental measurements of turbulence are crucial 
for us to understand it better, or at least its physical properties. For this rea- 
son, a 2D (8 radial x 4 vertical (poloidal) channels) beam emission spectroscopy 
(BES) diagnostic was installed on MAST to measure density fluctuations associ- 

^^Turbulence is intellectually so intriguing that proving existence and smoothness (no singu- 
larity) of the Navier-Stokes solution in three dimensional space is set to be one of the seven most 
important open mathematical problems by the Clay Mathematics Institute (www.claymath.org). 
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ated with turbulence. We know that the nonlinearity of a system and the linear 
drive of unstable modes are vital for the system to be in a turbulent state, other- 
wise the system is just filled with well-behaved (stable) waves. In this work, we 
restrict the description of plasma turbulence to its observed properties such as 
linear drives, spatial and temporal characteristics and nonlinearly saturated lev- 
els of turbulence (density fluctuations), and how these quantities are correlated 
with local equilibrium quantities. 



1.2.1 Drift waves and shear flows 

Tokamak plasmas are far away from a 'global' thermal equilibrium state, but 
rather in a force-balanced state: J x B = Vp with the plasma current density J 
and pressure p and the background magnetic field B. Note that plasmas are in 
thermal equilibrium 'locally'. Hence, there can exist gradients in the moments 
of phase-space distribution functions such as densities n, fluid velocities U and 
temperatures T. Gradients in density and temperatures of the ions Tj and elec- 
trons Te drive drift waves, which, in turn, can cause plasmas to be in a turbulent 
state via nonlinear self-interactions of drift waves Il9ll26]-l37|. More detailed de- 



scriptions of drift waves are provided in Appendix A 

While the gradients of density and temperatures are regarded as the linear 
drives, sheared mean plasma flows can both drive and suppress the turbulence: 
shear in the parallel component of the plasma flow Vf/y can drive turbulence 
p8] - |4T] l (as described in Appendix A.2[ >, while the perpendicular component of 



shear \/U± is known to suppress the turbulence ||42l - |52| . Here, f/y and U± are 
the parallel and perpendicular (with respect to the background magnetic field) 
components of the mean plasma flows, respectively. Thus, the ratio of the two 
shearing rates, Vf/y/Vf/j., plays an important role on the plasma turbulence in 
tokamaks [4l|. 

Let us digress briefly from the shear flow issue to introduce a new quantity 
called the 'safety factor' q. This will turn out to be one of the major parameters 
controlling the ratio Ww/VUx. The safety factor q is defined as BH 

where and are the toroidal and poloidal components of the magnetic field, 

respectively. The closed line integral {<f dl) is carried over one poloidal rotation at 

a fixed flux surface, q is interpreted as the number of required toroidal rotations 

for a magnetic field line to close itself, i.e., to complete one poloidal rotation. 

For a large aspect ratio tokamak with a circular poloidal cross-section, q can be 
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Figure 1.7: (a) A schematic of a helical magnetic field line with g = 4 on a flux 
surface. Figure from www.efda.org.; (b) A schematic of local magnetic field vec- 
tors (black arrows) in toroidal and poloidal Bq directions with the pitch angle 
a. Red arrows illustrate local mean plasma flows in parallel {U\\), perpendicular 
and toroidal {U^) directions. 



approximated as: 

where e = r/R (cf. = a/Ro). Figure LZf^a) shows an example of a flux surface 
with g = 4. Note that q is called the safety factor because of its close relation with 
the MHD (magnetohydrodynamic) stability: the higher the q, the more MHD 
stable it is BH. How q is related to the MHD stability is beyond the scope of this 
work. 

Suppose that the mean plasma flow is purely in the toroidal directiori[^ de- 
noted as U(p, then U\\/U± is equal to B^/Bq (see Figure [L7|[b)) which is, in turn, 
just q/e from Eq. (1.4 1. Thus, the ratio VU\\ /VU±, i.e., the ratio of turbulence drive 
and suppression, can be approximated as q/e, and we expect that the smaller q/e, 
the less turbulence. Experimental observations on the turbulence level as a function of 
q/e, one of the major topics of this work, is discussed in Chapter^ 



1.2.2 Spatial and temporal characteristics of plasma turbulence 

One of the common features in wide ranges of turbulent flows observed in 
nature, such as smoke from a cigarette, water flows in tunnels, oceanic currents, 
hurricanes, clouds, and the solar wind, is that they all have broadband spectra 
both in the frequency and the wavenumber domainj^ Plasma turbulence in a 

^■^The mean flows in the poloidal direction are strongly damped by the (neoclassical) collisions 
Il53ti56| , and Chapter|5]in this work shows experimental signatures consistent with the neoclassi- 
cal damping of the poloidal flow. 

^^However, the scalings can be quite different for different physical phenomena. 
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tokamak is no exception to this. Indeed, the plasma turbulence observed using 
the 2D BES diagnostic installed on MAST exhibits broadband frequency spectra 



(see Figure 3.5 ). Note that the BES turbulence data are not analysed in terms of 
wavenumber spectra due to the limited number of spatially separated channels 
of the diagnostic!^ The principle of this diagnostic, without which this work would not 
exist, is explained in Chapter^ 

The broadband spectra of turbulence means that there exist many differ- 
ent spatial and temporal scales from the energy-containing (or energy-injection) 
range down to the dissipation range through the inertial range if these ranges are 
well separated in their scales. In this work, we concentrate on the spatial and tem- 
poral scales of the plasma turbulence in the energy-containing range, k^pi < 1, 
due to the limited spatial resolution of the 2D BES diagnostic, i.e., the sensitiv- 
ity of the diagnostic decreases for the higher wavenumber ranges. Here, k± is 
the perpendicular wavenumber of the turbulence and pi the ion Larmor radius. 
Nevertheless, because the energy-containing range gives the largest contribution 
to the turbulent transport due to the large scale of the fluctuations, studying and 
understanding the turbulence characteristics in this range is the most critical to 
building a more efficient fusion power plant. 

Spatial structures of the plasma turbulence in a tokamak can be thought about 
in terms of time scales, a consequence of the v ■ V operator which arises in the 
evolution equations where v can be related to various physical effects. For in- 
stance, f I V|| acting on some turbulence distribution function h in phase-space is 
associated with the parallel dynamics of the turbulence, hence the parallel struc- 
ture of the turbulence. Thus, it can be identified what physical effects influence 
(or possibly determine) the spatial structures of the turbulence by comparing var- 
ious time scales in a tokamak. In fact, we do find from our experimental data that 
the parallel streaming time, the magnetic drift time and the linear drive time associated 
with drift waves driven by density- or ion-temperature-gradient scale with the turbu- 
lence correlation time consistently. These results together with the consequences of these 
"balanced time scales" on the turbulence spatial structures are discussed in Chapter^ 



^*We may be able to construct wavenumber spectra from the data using the maximum entropy 
method Il57ll58l , but such a technique has not been used. This is left as a future work. 
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Chapter 2 



Structure of this work 



Our goal is to observe and understand the turbulence in MAST, hence we 
installed a 2D BES diagnostic on MAST, and this work is based on the experi- 
mental data from the diagnostic. We structure our work in the following logical 
steps: 

(1) We need to obtain useful information, such as the properties of turbulence, 
from raw BES data. Thus, we develop statistical techniques to do so in Chapter 

13 

(2) The statistical techniques must be examined for their reliability, hence we 
develop a way to generate synthetic 2D BES data in Chapter |4] for which all the 
properties of turbulence are known in advance so that the statistical techniques 
can be examined. 

(3) Then, we obtain turbulence information from the 2D BES data using the 
statistical techniques, and make conclusions about physics of turbulence in 
subsequent chapters: what causes fluctuating density patterns to move in the 
poloidal direction in the lab frame while the plasmas are rotating toroidally 
(Chapter |5]), and what equilibrium quantities are correlated (or determine) 
spatial and temporal characteristics of turbulence (Chapter [6]). 

(4) Finally, we describe how we can improve the performance of tokamaks with a 
careful correlation analysis among many local equilibrium parameters in Chap- 
ter [7| by producing a critical manifold separating turbulent and non-turbulent 
state. 

We first describe in more detail how the 2D BES diagnostic can be used 
to measure density fluctuations up to a few 100 kHz range with the spatial 
resolution of ~ 2 cm, and how we obtain the characteristics of turbulence from 
the BES data in Chapter |3j Then, we explain in Chapter |4] how to generate 
synthetic BES data, a forward model of the 2D BES diagnostic, which can be 



used in many different aspects such as examining the reliability of statistical 
analyses as in Chapter |5] and performing direct comparisons between outputs of 
numerical turbulence simulations and experimental turbulence data Il59i - |6T1 . 

We find that the fluctuating density patterns move in the poloidal (vertical) 
direction in the lab frame. This seemingly contradicting result with the mean 
plasma flows being dominantly toroidal is discussed in Chapter |5] We show, 
via a careful ordering of the density continuity equation, that such an apparent 
poloidal motion of the patterns arises due to the fact that the elongated patterns 
in the parallel direction are advected by the dominant mean toroidal plasma 
flows provided there exists a finite angle between the parallel and toroidal 
directions, i.e., the projection effect analogous to the apparent up-down motion 
of helical strips of a 'rotating barber-pole' [|62|. 

Next, we compare the turbulence correlation time, the particle (ion) parallel 
streaming time, the drift time associated with ion temperature or density gradi- 
ents and the magnetic drift time finding that they are all comparable in Chapter 
|6j This result suggests that the observed turbulence in MAST is "critically 
balanced" and its characteristics are determined by the local equilibrium, from 
which spatial correlation lengths of turbulence are derived and examined. Fur- 
thermore, we infer the turbulence nonlinear time from the density fluctuations, 
and we find that the ratio of the inferred nonlinear time to the turbulence corre- 
lation time is a function of ion-ion coUisionality We argue that this observation 
is consistent with the decorrelation of turbulence being dominantly controlled 
by zonal flows H63|l64l. 

The final result we present in this work provides a way to achieve better 
tokamak performance. We show in Chapter |7] statistically that the normalized 
ion-temperature-gradient scale length R/Lt^ is inversely correlated with local 
q/e, the ratio between the shearing rates of parallel and perpendicular flows; 
while R/Lxi and the local shearing rate of the mean toroidal plasma flows are 
positively correlated. The dependence of R/Lt^ on q/e is strong which implies 
that R I Lx^ can be increased at a fixed shearing rate of the mean flows by lowering 
q/e. Furthermore, we present a critical manifold in local equilibrium parameter 
space separating turbulent and non-turbulent state based on the fact that the 
observed turbulent heat flux is inversely correlated with R/ Lt^- 

Then, we close this work with the conclusions in Chapter |8j 

The author of this work generated point-spread-functions of the 2D BES 

system, developed software to generate synthetic 2D BES data and to perform 

statistical analyses on raw BES data. A. Field, D. Dunai and S. Zoletnik de- 
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signed, installed the 2D BES system on MAST and took the data. C. Michael 
generated EFIT data which contain the equilibrium magnetic field information. 
A. Schekochihin, I. Abel, M. Barnes, G. Colyer, S. Cowley, F. Parry, E. Highcock 
provided theoretical insights. 
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Part II 

Turbulence (Density Fluctuation) 
Measurements 



Chapter 3 

Beam Emission Spectroscopy (BES) 
diagnostic: Measuring density 

fluctuations 



This chapter is largely based on Refs. I[65\\66t 

3.1 Principle of 2D BES turbulence diagnostic 

The 2D BES system on MAST utilizes an avalanche photodiode (APD) 2D ar- 
ray camera [i67il with eight columns and four rows of channels, which have an 
active area of 1.6 x 1.6 mm^ each. It measures the Doppler-shifted emission 
from the collisionally excited neutral-beam atom^ (deuterium) with a temporal 
resolution of 0.5 /is. The optical system is designed such that the observed loca- 
tions can be radially scanned along the path of the neutral beam (South-Neutral 
Beam Injection (S-NBI)) whose 1/e half -width is 8cm, while the optical focal point 



follows the axis of the beam (see Figure 3.1 1. More detailed descriptions on the 



optical system of the 2D BES are available in Ref. [i68i . The nominal location of 
the BES system, i.e., where the optical line-of-sight (LoS) is best aligned with the 
local magnetic fieldj^is at major radius R = 1.2 m. At this location, a magnifi- 
cation factor of 8.7 at the axis of the beam results in each channel observing an 
area of 1.5 x 1.5 cm^ with 2 cm separation between the centres of adjacent chan- 



nels (see Figure 3.2). The poloidal (vertical) locations of the views are fixed at 



^Neutral beams are injected mainly to heat and to provide toroidal momentum to MAST 
plasmas. 

^The BES measurements are spatially localized to a good degree because the LoS at the inter- 
section of the beam are approximately parallel to the local magnetic field lines and the intersection 
length is much shorter than the parallel correlation length of the turbulence. 
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R 




S-NBI 



-2 m 2m ^ 

Figure 3.1: Top view of MAST. 
South-NBI is injected to heat and 
rotate (toroidally) the plasmas (red 
strip), and the optical system of 
the 2D BES is focused to the axis 
of the S-NBI (blue rectangle) to 
collect Doppler-shifted emis- 
sions from the collisionally excited 
neutral-beam atoms. The 2D BES 
can be scanned radially following 
the axis of the S-NBI. A typical in- 
jection energy of the S-NBI is 60 — 
70 keV with the 2 - 3 MW of input 
power. 




mm) 



Figure 3.2: Side view of MAST. Ra- 
dial viewing position of the 2D BES 
can be scanned from the core to the 
edge of plasmas. The 2D BES sys- 
tem consists of 8 radial x 4 poloidal 
channels with 2 cm separation be- 
tween the neighboring channels in 
both directions. The area of each 
channel is about 1.5 x 1.5 cm^ which 
comes from the magnification factor 
of 8.7. This varies slightly depend- 
ing on the radial viewing position of 
the 2D BES system. 



Z = -0.03, -0.01, 0.01 and 0.03 m. 

The angle between the LoS of the 2D BES system and the neutral beam with 
the injection energy of 60 — 70 keV results in a Doppler shift of the Da emission 
approximately 3 nm to the red from the background Da at 656.28 nm (see Fig- 
ure 



'5Mn The background Da can be removed with a suitable optical filter [68J, 



and so that the Da emission detected by the 2D BES system comes only from the 
neutral beam, hence the measurement is localized to the beam. A more detailed 
description of the 2D BES system on MAST can be found elsewhere |I68|[. 

The measured intensity of the Da beam emission is directly related to the 
background plasma density because the latter is the cause of the excitation of the 
neutral-beam atoms. The beam atoms are excited by collisions with the electrons, 
ions and impurities, but at energies greater than 40 keV, the electron contribution 
can be ignored ||69|. The fluctuating part of the plasma (ion) density 5n can be 

^There exists strong background Da emission from the plasma edge. 
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determined according to 

n Pbes I 

where n is the mean plasma density, and SI and / denote the fluctuating and 
mean parts of the photon intensity, respectively. /3bes is a coefficient depending 
on the population of the excited state and is a weak function of the background 
plasma density with values in the range 0.3 < /3bes < 0.7 with the negligible tem- 
perature dependence ||70|. /3bes is calculated according to a coUisional-radiative 
model, and we use values from the work of Hutchinson in [701. A careful mod- 
eling of the 2D BES system shows that the system is capable of measuring the 
density fluctuation levels down to a few 0.1 % [68] 

Thus, the 2D BES system on MAST directly measures fluctuations of plasma 
density in the radial (8 channels)-poloidal (4 channels) plane at a fixed toroidal 
location with the spatial resolution of ~ 2 cm in both directions when the S-NBI 
is injected]^ Note that the actual spatial resolution is somewhat broadened due 



to various physical effects which are explained in more detail in Section 3.2 



Figure 3.3 shows the evolution of (a) the plasma current, (b) line-integrated 
electron density, (c) S-beam voltage and (d) BES intensity data from three differ- 
ent radial locations, R = 123 cm, 119 cm and 115 cm, at a fixed poloidal location 
Z = 1.0 cm. It shows that the levels of BES signals rapidly increase after the S- 
NBI is injected at t = 0.05 s as the size of plasmas increases until around 0.1s. 
This increase in the size of plasmas can be seen from Figure 3.3 'e)-(f) showing 
the contour of normalized flux surfaces ip, i.e., = 1.0 at the LCFS, att = 0.05 s, 
0.075 s and 0.1 s. The circles show the viewing locations of the BES system (three 
green circles correspond to the positions where the BES signals are obtained for 
Figure 3.3[ d)). Then, the BES signals slowly increase as the electron density in- 
creases until the S-NBI cuts off at t = 0.2 s. Figure 3.4 shows measured spectra 
of beam emission from Z = cm at the same toroidal location as the BES system 
using the Motional Stark Effect (MSE) diagnostic viewing at a radial location of 
R = 1.2 m. The observed Doppler shifts of the Da lines are almost the same as 
those for the BES system. Spectra are shown during a period with the S-beam 
voltage on (red) and during the ohmic phase just after the beam is switched off 
(black) (see the yellow strips in Figure |3.3[ d)). There are three Doppler shifted 
Da lines corresponding to the emissions from D, D2 and D3 corresponding the 
full, half and third energy components of the neutral beam. There are non-zero 



^Typical fluctuation levels at the edge of plasmas are of the order of 1 % and 0.1 % in the core 
for L-mode discharges, while core regions in H-mode discharges usually have smaller fluctuation 
levels. 

^BES signals are not spatially localized when S-NBI is not injected (see Figure 3.3 d)). 

18 



3.1. Principle of 2D BES turbulence diagnostic 



Flux surfaces 



1.0 



(a) Plasma current 



_ 1.5 

cjj ' 

1.0 

O 

" 0.5 
0.0 



(b) Line-integrated density 



(c) S-beam voltage 




R=123cm 
119cm 
1 15cm 




20 



E 
u, 

N 




-20 



20 



E 
o, 

N 



-20 



20 



E 
u, 

N 



-20 




, 03 / ' CV 



0.2 0.3 
Time [sec] 



0.4 



0.5 



80 100 120 140 
R [cm] 



Figure 3.3: Evolution of (a) the plasma current (b) line-integrated electron density 
(c) S-beam voltage and (d) BES intensity data at R = 123cm (Ch.#19: black), 119cm 
(Ch.#21: red) and 115 cm (Ch.#23: blue). The vertical green dashed lines show the 
times when the 2D BES system obtains localized density fluctuations from the 
S-NBL The two yellow strips (0.18 - 0.20 s and 0.22 - 0.24 s) indicate the time 



durations when the spectra of beam emission are measured shown in Figure 3.4 
(e), (f) and (g) show the contour of flux surfaces ijj att = 0.05 s, 0.075 s and 0.1 s, 
respectively. Circles show the viewing positions of the BES system (the green 
circles indicate the locations where BES signals are obtained for (d)). 



levels of BES signals even after the S-NBI cut-off (Figure 3.3[ d)) which is due 
to the CII (657.81, 658.29 nm) lines as the carbon impurity density is not negligi- 
ble in MAST. Figure |3.5| shows a typical cross-power spectrum and cross-phase 
between the two poloidally separated channels with and without MHD modes. 
Effects of MHD modes in determining turbulence characteristics are discussed in 
Section l33l 
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Figure 3.4: Spectra of beam emission measured both before (red) and after (black) 
the S-beam cut-off time using the Motional Stark Effect (MSE) spectrometer at a 
viewing radius of 1.2 m together with the transmission of the BES filter (blue). 
Figure taken from Ref. ||68|. 

3.2 Point-spread-functions of the 2D BES system 

BES is a volume-sampling diagnostic, therefore in order to generate synthetic 
BES data (see Chapter |4]) one really requires LoS integration based on three- 
dimensional input data of the fluctuating plasma density which would require 
very large data files. However, by availing the spatial structure of the turbu- 
lence which is elongated along the magnetic field lines synthetic BES data can be 
constructed with 2D input data on the poloidal cross-section using the 2D point- 
spread-functions (PSFs) of the detectors. The relevant physical effects such as 
magnetic field-line curvature, LoS geometry, finite excited-state lifetime of Da, 
the beam attenuation and divergence must all be taken into account in the cal- 
culation of the PSFs. Without considering these effects, the perpendicular spatial 
resolution will be limited to ~ 2 cm. 

In order to calculate the 2D PSFs, successive image planes are constructed 
along the LoS |7l[. The size of the detector images at these planes is set by the 
optical magnification factors which vary along the LoS. Light-cones whose sizes 



20 



3.2. Point-spread-functions of the 2D BES system 




100 200 300 400 
Frequency [kHz] 



500 




100 200 300 400 500 
Frequency [kHz] 



Figure 3.5: Spectra of crosspower (top panel) and crossphase (bottom) of a mea- 
sured fluctuating density from 2D BES system on MAST when the signal contains 
MHD modes (red dashed) or no MHD modes (black solid). 



are determined by the optical system are convolved with these detector images. 
The resulting blurred images are then convolved with an exponentially decaying 
function whose 1 /e distance is calculated by considering the beam velocity pro- 
jection of the perpendicular to the LoS and the half -life of Da- Hutchinson |[70| 
calculated the half-life of denoted as r^, and it is reported that at a plasma 
density of ~ 10^^ m~^, is ~ 3 — 10 ns. The resultant images along the LoS are 
then moved to the optical focal plane by following the magnetic field-lines. 

Figure 3.6 shows examples of generated PSFs at two different radial locations 
covering (a) 1.00 m < R < 1.14 m and (b) 1.16 m < R < 1.30 m. The radial smear- 
ing effects are mainly due to the finite half -life of Da as the excited neutral-beam 
atoms can travel (radially) finite distances before they emit Da fluorescence. On 
the other hand, the poloidal smearing is due to the pitch angle a of the magnetic 
field. The pitch angle changes much more significantly from the core to the edge 
of plasmas in a spherical tokamak compared to a conventional tokamak, and this 
effect can be seen in the strong dependence of the poloidal width of the PSFs as a 
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Figure 3.6: Example of PSFs of the 2D BES system covering (a) 1.00 m < R < 
1.14 m and (b) 1.16 m < R < 1.30 m. The lines (red and blue) are the 1/e contour 
lines of the PSFs, and the asterisks show the optical focal points. 



function of viewing radius. Consequently, deconvolution of the measured signal 
using the calculated PSFs on MAST is not trivial as the functions vary in space 
which is not the case in a conventional tokamak such as in DIII-D tokamak ||7T|. 
More sophisticated algorithms will be required to enable the deconvolution with 
spatially varying PSFs [[72l|. Utilizing such an algorithm is beyond the scope of 
this work. In fact, we find that most of measured turbulence has poloidal corre- 
lation lengths longer than the poloidal width of PSFs, and we do not use the data 
points whose radial correlation lengths are less than 2.0 cm. 
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3.3 Extracting turbulence characteristics from the 
BES data 

Once the raw BES data are obtained, we need to extract useful information, 
i.e., spatial and temporal characteristics of turbulence, from the data. This is done 
by calculating the spatio-temporal covariance function and its normalized ver- 
sion, the correlation function, throughout this work. The covariance function is 
calculated as 

C {AR, AZ, At) = {61 {R, Z, t) 6I{{R + AR,Z + AZ, t + At)) , (3.2) 

and the correlation function 

^ / A ^ A ^ A X (AR, AZ, At) 

C (AR, AZ, At) = ^ ' ' ^ (3.3) 

y/{SP (R, Z, t)) {SP {R + AR,Z + AZ, t + At)) 

where R and Z denote the radial and poloidal (vertical) coordinates with their 
channel separation distances A_R and AZ, respectively, and t the time and At the 
time lag; and (■) denotes a time average. Note that we use time-averaging in- 
stead of ensemble-averaging throughout the work by assuming that our signal is 
ergodicj^ 

One must be careful obtaining the fluctuation level of the BES signal as 
the information is acquired from an auto-covariance function at At = 0, i.e., 
(5/2 (t)) = (Ai? = 0,AZ = 0, At = 0), which contains the noise of the BES sig- 
nal as well. To remove the noise from the signal, we independently measure pho- 
ton noise levels by illuminating the BES sensors with an LED calibration source]^ 
We obtain 150 different DC levels of noise signals from to 1.5 V creating a 
database of noise signals. Let us define a measured BES signal / (t) from a MAST 
discharge consisting of a mean (/), plasma turbulence signal S (t) and noise (t): 
I (t) = (/) + 51 (t) = {I) + S (t) + N{t). The true (squared) fluctuation level of the 
turbulence is {S^ (t)) = {AR = 0,AZ = 0, At = 0)~{N^ {t)) assuming that S {t) 
and {t) are uncorrelated. As the noise source of the 2D BES system on MAST 
is dominated by the photon noise ||67L we find the corresponding {t) from the 
created noise database with the LED light whose DC level matches with the DC 



level of a measured BES signal, i.e., (/). This is shown in Figure 3.7 'a) where the 



^Ideally, we want to repeat the same experiment many times and average the data from the re- 
peated experiments. However, as repeating the experiment is expensive (not to mention whether 
repeating it is possible or not), such ensemble-averaging is not possible in practice. 

''An LED calibration source is mounted inside the optics box of the BES system, so that we 
can apply LED light to the BES channels while all the other light sources are blocked by closing 
the shutter of the optics box. 
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Figure 3.7: (a) Time trace of measured raw BES signal (black) and measured pho- 
ton noise (red) whose mean voltage matches with that of the BES signal; (b) Cal- 
culated auto-covariance functions of the BES signal (black), the photon noise sig- 
nal (red) and the photon noise subtracted BES signal (blue). Note that both BES 
and photon noise signals are filtered from 20 to 100 kHZ in (b). 



black line is a measured BES signal and red dashed line is the corresponding pho- 
ton noise with the same DC levels. If one wants to frequency-filter the measured 
signal, then it is necessary to filter the (t) with the same frequency band as 
well. Figure |37[b) shows the C {AR = 0,AZ = 0, At) (black), {N (t) N{t + At)) 
(red) and C {AR = 0,AZ = 0, At) - {N (t) N {t + At)) (blue) where both SI (t) 
and (t) are frequency-filtered from 20 to 100 kHz. The fluctuation level of S (t) 
is, then, the square-root of the {AR = 0,AZ = 0, At = 0) - {N^ (t)), i.e., the 
square-root of the blue curve at At = 0. Note that we remove beam noise by 
frequency filtering the signal below 20 kHz. 

Spatial correlation length is estimated using the correlation values at At = 
with finite values of AZ or AR. We fit C (AR = 0, AZ, At = 0) to the function 
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Figure 3.8: Examples of C {AR = 0, AZ, At = 0) from MAST shot #27267. Evo- 
lution of (a) toroidal plasma flow and (b) cross-power spectrogram of 
BES signal between two poloidally separated channels by 2 cm; (c) shows 
C {AR = 0, AZ, At = 0) from the time interval (1) in black and (2) in red. Time 
intervals (1) and (2) are marked in (b). Note that the time interval (2) contains 
MHD activities (starting at t ~ 0.23 s with multiple harmonics) which forces pz 



in Eq. (3.4) to be large from the fitting procedure. In both cases, poloidal correla- 
tion lengths £z are estimated to be 10.0 cm. 



fz {AZ) defined as 

fz (AZ) =pz + {l- Pz) cos 



2tt- 



AZ 



exp 



\AZ\ 



(3.4) 



to estimate the poloidal correlation length iz- Here, we assume wave-like fluc- 
tuations in the poloidal direction ||73| with the same wavelength and correla- 
tion length. It is not possible to distinguish between the two with only four 
poloidal channels. However, this assumption is verified using the 2D BES data 
from DIII-D tokamak which has 8 poloidal channels]^ The constant pz is a fit- 
ting parameter to account for global structures such as MHD modes. This is 
necessary as 2D BES data from MAST discharges usually include MHD sig- 



nals (see Figure 3.8 'b) and Section 5.4.1 1. Consider a fluctuating part of BES 
data 6Ii (t) = Si (t) + Gi (t) + Ni (t) where S, G and denote for turbulence, 
global and noise signals, respectively. The subscript i stands for a BES channel 
number. Assuming that these three signals are not correlated to each other and 
noise signals from different channels are uncorrelated, we have {5Ii (t) 51 j (t)) = 
{Si (t) Sj (t)) + {Gi (t) Gj (t)). With a fixed subscript i and varying j in the poloidal 

^The author of this work visited General Atomics where DIII-D tokamak is located in May 
2012 and verified the assumption with Dr. George McKee who is in charge of its 2D BES system. 
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direction, (Gj (t) Gj (t)) is a constant, i.e., pz in Eq. (3.4 1, if there exists a global 
mode whose spatial structure is larger than the poloidal extent of the BES sys- 
tem; while {Si (t) Sj (t)) is assumed to have a wave-like spatial structure. Figure 



3.8 shows examples of the fitting results: one from an MHD quiet period marked 



as (1) in (b) and the other from an MHD active period marked as (2) in (b). The 
fitting gives pz = 0.04 and pz = 0.90 for these two cases, respectively, as shown 
in Figure 3.8 'c). The estimated £z is 10.0 cm in both cases. As a large pz is likely to 
affect the estimation of iz, and uncertainty levels are larger when MHD modes 
are active, we set data points with large pz unreliable and do not use them for 
the purpose of turbulence studyj^ One of the many criteria we use to reject data 
points with MHD modes is pz > O.sj^ Finally, we note that the BES cross-power 
spectrogram shows a drop of turbulence level when plasmas start to rotate sig- 
nificantly in the toroidal direction shown in Figure [3^ a) and (b). The influence 
of such rotation on the turbulence level is a subject of Chapter [7| 

The radial correlation length is estimated by fitting C ( Ai?, AZ = 0, At = 0) 
to the function (AR) defined as 

r I Ai?i 

fn (AR) =pji + {l- pr) exp 



(3.5) 



where pr plays the same role as pz did for fz- Note that we do not use a wave- 
like structure in the radial direction as observed in other tokamaks Il60l[73ll . This 
is also vindicated for measured turbulence in MAST as C (AR, AZ = 0, At = 0) 
usually monotonically decreases from AR = 0, and their values are non-negative 
as in Figure 3.9 Statistically, pr and pz should be more or less the same, but they 
can differ as estimating a constant offset in /r (AR) is non-trivial, i.e., there is no 
obvious baseline. Thus, we estimate Ir only for those points with pz < 0.5. If the 
fitting results have pr > 0.5 even for the points with pz < 0.5, then we cut-out 
these points as well. In other words, if a data point satisfies pz > 0.5 or pr > 0.5, 
then we set the point unreliable. Figure 3.9 shows the results for the two cases 
mentioned earlier: without MHD modes in black and with MHD modes in red. 
For both cases, pr is estimated to be 0.0 with Er = 6.9cm without MHD and 2.5cm 
with MHD. Again, it is clear that uncertainty levels are worse when there exist 
MHD modes. 

In addition, a spurious long range radial correlation can be introduced to the 
signal due to edge-induced fluctuations known as beam common-modes [74J. 
The beam common-mode can be identified by looking at the phase difference be- 
tween the signals obtained from the core and edge of plasmas, and we find little 

'However, these points are not only valid but also valuable for the study of MHD activities. 
^"Other criteria are listed in Chapter |6j 
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Figure 3.9: Examples of C (Ai?, AZ = 0, At = 0) from the same time intervals as 



in Figure 3.8 The fitting procedure gives the radial correlation lengths Er of 6.9cm 



and 2.5 cm for interval (1) and (2), respectively. When MHD activities are present, 
the fitting is much worse than that without any MHD activities. 



signatures of such modes in our BES data. 

Estimating the correlation time Tc of measured turbulence has to take 
account of the fact that turbulence is advected in the lab frame due to 
toroidal plasma flows. Thus, we use the Lagrangian approach ||74|. We 

fit C {AR = 0, AZ, At = T^l^^ (AZ)), where t^I^^^ (AZ) is the time delay when 
C {AR = 0, AZ, At) is maximum at a given AZ, to the function (AZ) defined 
as 



fr (AZ) = exp 



peak 



Tr 



(3.6) 



Figure 3.10 shows examples of the fitting results: (a) without MHD and (b) with 
MHD modes. The cyan lines show the fr (AZ) fit with (a) Tc = 11.8 /xs and 
(b) Tc = 2.4 fxs. Again, when MHD modes are present, uncertainties are larger, 
and it can potentially affect the estimation of Tc. The described Lagrangian 
approach relies on the fact that how precisely r^l^-^ i'^Z) can be located; how- 
ever, the existence of MHD modes can modify the true Tp^g^,^ i^Z) significantly. 



an effect which is described in Section 5.3 in more detail. Finally, we note 
that the reliability of this method depends on the temporal decorrelation dom- 
inating over the parallel spatial decorrelation because we would see decreasing 
C {AR = 0, AZ, At = Tp^g^j^ i^Z)) as a function of AZ even if Tc approaches infi- 
nite if the parallel correlation length is not long enough. This fact is quantitatively 
discussed in Chapter |6] 
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Figure 3.10: Examples of C (Ai? = 0, AZ, At) from the same time interval as in 



Figure |3.8l Circles indicate the positions of t^I^^ (A^) and cyan lines are the 
fr (AZ)W with (a) = 11.8 ^is and (b) = 2.4 /is. Note that t^^^ (AZ) is heavily 
influenced by the MHD activities in (b) which causes (possibly) an ambiguous 
estimation of 
in Section l53l 



The effect of MHD modes on ( A2^) is studied in more detail 



All the fits described in this section are obtained via the mpf it procedure [75J. 
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Chapter 4 
Generating synthetic BES data 



This chapter is largely taken from Ref. Ii66t 

As discussed in Chapter |3| deconvolution of the measured signal using the cal- 
culated PSFs are not trivial. However, we can generate synthetic BES data using 
the PSFs, which then can be used to find any effects of PSFs on statistical analyses 
as in Chapter [s] and to make direct comparison^ of numerically simulated data 
with the experimental data as in Ref. HSU. 

Note that we use a local coordinate system with x denoting radial direction 
and y denoting poloidal direction in this chapter. 



4.1 Gaussian eddies in space and time 

To be able to numerically examine the reliability of any statistical analyses on 
turbulence data, it is necessary to know the exact characteristics of the turbulence. 
For this purpose, we numerically generate artificial fluctuating density patterns, 
random both in space and time, then produce synthetic BES data using the PSFs 



of the 2D BES system on MAST (as described in Section 4.2 1 and compare the 
inferred characteristics from statistical analyses with the true characteristics. 

We follow a similar approach to the one suggested by Zoletnik et al. [76J. Let 
the density patterns be described by Gaussian structures both in space and time. 



^ These results are not included in this work partly because the work is mainly carried by Dr. 
Anthony Field using the calculated PSFs described in Chapter |3] 



4.1. Gaussian eddies in space and time 



namely. 



N 

5n (x, y, t) = 2, '^^oi exp 




1=1 



cos 



(4. 



where x, y and t denote radial, poloidal and time coordinates, respectively. These 
numerically generated density patterns are referred to as "eddies" in this work. 
Here is the total number of eddies and the subscript i denotes the i^^ eddy 
in the simulation; SriQi, xoi, yoi and toi are the maximum amplitude and central 
locations in the x, y and t coordinates of the i*^ eddy, respectively; X^r \ and 
Tiife are the widths of our Gaussian eddies in the x, y and t directions; Tufe is the 
lifetime (or the correlation time) of the eddies in the moving frame; Vy{t) is the 
apparent advection velocity of the eddies in the poloidal direction. Although it 
is possible to introduce a finite radial velocity shear by making Vy a function of x, 
the effect of such shearing rates is not investigated in this work, so we will only 
consider Vy that are independent of x. The cos term in the y (poloidal) direction 
is introduced to model wave-like-structured eddies in the poloidal direction as 
observed in tokamaks II73I1 . In choosing the wavelength in the poloidal direction, 
we assume that it is the same as the exponential decay length]^ This is because the 
two lengths cannot be separately and reliably measured using only four poloidal 
channels. Note that the envelope (i.e., the exp term) and the wave structure (i.e., 
the cos term) of 5n{x, y, t) have the same advection velocity Vy{t). The central lo- 
cations of eddies, xq, yo and to, are selected from uniformly distributed random 
numbers, whereas their amplitudes Sno are selected from normally distributed 
random numbers whose standard deviation is oneE] 

The spatial domain of the simulation is 25 cm and 20 cm with the mesh size of 
0.5 cm in radial (x) and poloidal (y) directions, respectively. The time duration of 
the simulation is 20 ms with a 0.5 /is time step so as to have the same Nyquist fre- 
quency as the real 2D BES data from MAST. The widths A^. and A^, are set so that 
the full width at half maximum (FWHM) in the radial direction and the wave- 
length in the poloidal direction are ~ 8 cm (i.e., A^ = 3.53 cm) and ~ 20 cm (i.e., 

^The author of the work visited General Atomics and confirmed that this assumption is rea- 
sonable with the 8 channel poloidal measurements on DIII-D. 

■'it is worth mentioning that there is another scheme of generating such eddies numerically, 
proposed by Jakubowski et ah llTTl . They generated the time series of fluctuating density (Sni) 
using the inverse Fourier transform of a broadband Gaussian amplitude distribution in frequency 
space. Then, a second signal {6n2) was generated by imposing the desired time-delay fluctuation 
on the Sni such that 6n2 was a time-delayed version of Sni. This method does not include spatial 
information for the signals. 
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Xy = 20.0 cm), respectively, which are similar to the measured correlation lengths 
with the 2D BES system on MAST|^The eddy lifetime in the moving frame {rufe) 
is set to 15 us. However, some of the data sets in this work have different val- 
ues of Tii fe, so the effect of tu on a statistical analysis, the cross-correlation time 
delay (CCTD) method to measure the velocity of fluctuating patterns, can be in- 



vestigated in Section 5.3 



The total number of eddies is = 20000. If the eddies are too sparse in 
the simulation domain, then we may not achieve steady statistical results, while 
overly dense eddies may cause an effective widening of the specified spatial (A^ 
and Xy) and temporal {rufg) correlations as many eddies can merge into one larger 
eddy. Thus, we introduce another control parameter, the spatio-temporal filling 
factor (F), defined as 

XxXy \ f T[i 



N ■ ■ — , (4.2) 

total simulation area / V total simulation time/ ' 



where r^c is the autocorrelation time calculated as 1179 II 

(4.3) 



riife (K/Vy) 



• ac I 

y'Tufe + iK/^y 

for the generated eddies defined by Eq. ( 4T] >. All of our synthetic data was gen- 
erated so as to F ~ 0{1). 

The testing of the CCTD method in Section |53 will involve exploiting what 
happens if Vy{t) has a mean and a temporally varying components. Thus, we 
generate a temporal structure of Vyi at each x, 

Vy it) = {Vy) + SVy (t) 



(Vy) + Vy (t) * eXp 



'flue 



sin {27^ ffiuct) (4.4) 



where {vy) and 5vy are the mean and temporally varying velocities, respectively, 
Tfiue and ffiuc are the lifetime and frequency of 6vy{t), respectively, and Vy{t) is 
generated from normally distributed random numbers. The RMS (root-mean- 
square) value of 6vy{t) denoted as Sv^^^ will be varied as well as {vy) to inves- 
tigate the effects of these quantities on the CCTD method. r/;„c and ffiue allow 
one to introduce structured temporally varying velocities, while the randomness 
is kept by Vy. As one of the causes for the temporal variation of the poloidal ve- 
locity is believed to be the existence of geodesic acoustic modes (GAMs)|^[80J, we 

^Note that Smith et al. ITS! also reported that poloidal correlation lengths of the density pat- 
terns are ~ 20 cm using their 2D BES system on NSTX. 

^We do not investigate whether the CCTD method is able to detect such a temporally struc- 
tured 5vy [t) (or GAMs) in this work, rather we investigate how the existence of these structures 
affects the CCTD-determined mean velocity. 
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choose Tfiuc = 500 and ffi^c = 10 kHz to mimic the GAM features detected by 
Langmuir probes on MAST |8l|. 

The simulations have been run on a NVIDIA(r) GeForce GTS 250 GPU card 
using CUDA programming, which increases the computational speed owing to 
the highly parallelizable structure of Eq. (4.1 1. 



4.2 Synthetic 2D BES data 

We generate the i*'* (1 to 8) radial and j*'' (1 to 4) poloidal channel of the syn- 
thetic BES data (t) by using the calculated point-spread-functions (PSFs) of the 
actual 2D BES system on MAST 1165 II described in Chapter |3] and Sn (x, y, t) from 



Eq. (4.1 1 with an additional random noise. In general, 6n (x, y, t) can be taken 
from any turbulence numerical simulations as done in Ref. ||59]1. Furthermore, a 
large-scale (in space) coherent (in time) oscillation is included to imitate a global 
MHD mode. Namely, (t) is defined as 

(t) = Itc + ^I'' it) + itiHD it) + 1% it) , (4.5) 

where ri)(j is the DC value - a typical value of 0.8 V is used for all channels [68J . 
The rest of the terms are as follows. 

6P^ (t) is the fluctuating part of the signal generated from the fluctuating den- 
sity 5n (x, y, t), i.e., the Gaussian eddies given by Eq. ( |4.1[ > in this work, convolved 
with the PSFs of the 2D BES system: 

./«W^«--'//^n(..,.OP«(..,)d.d,. (4.6) 

where (x, y) is the PSF of the i*^ and channel of the 2D BES system, nor- 
malized so that RMS value of SP^ (t) is 51^^^^. This value is set so that the ratio 
of 5J«*'^^ to I'ic is 0.05. 

Imhd (^) iTiodels an MHD (global) mode. We assume that the spatial scale 
of the MHD modes is larger than the BES domain in the poloidal direction, so 
^MHD i^) does not vary in the poloidal direction. The model MHD signal is gen- 



erated in a way similar to temporal behaviour of Vy (t) using Eq. (4.4), except that 
the mean value of II^j^d (^) zero and Tfi^c = 250 /xs. The frequency of the mode 
fAiHD and its RMS value, denoted Im^^, will be varied in various tests. The 
value of Tfiuc here is representative of MHD burst-like fishbone instabilities [|82| 
or chirping modes ||83ll in tokamaks, for which the spectrum has a finite band- 
width. 

Pj^ (t) represents the noise in the signal. As the noise of the 2D BES system on 

32 



4.2. Synthetic 2D BES data 




Figure 4.1: Autopower spectra of synthetic 2D BES data for various {vy). Note 
that the spectrum for {vy) = 10.0 km/s (green dash dot line) has finite I^mhd (i.e., 
temporal oscillations due to global modes) in Eq. (4.5 1 at 15 kHz, with fluctuation 
level of 5 % of the DC level. For other cases, If^o = 0. All the spectra are 
generated using a high-pass filter with the frequency cutoff at 5 kHz. 



MAST is dominated by the photon noise II67L (t) is generated using normally 
distributed random numbers. Its RMS level is set such that the signal-to-noise 
ratio (SNR) is 300, which is typical of the 2D BES system on MAST Il68ll . 

Figure |4^ shows examples of autopower spectra of the synthetic 2D BES data 
for {vy) = 2.0, 5.0, 10.0 and 40.0 km/s. The autopower spectrum is calculated as 
\FT {P^ (t)}\ where FT {■} is the Fourier transform in the time domain. Increas- 
ing the value of {vy) has two effects: Doppler shift and broadening of the spectra, 
as expected. Note that in Figure 4.1, the data for (vy) = 10.0 km/s contains the 
finite Jf/^f, with Jmnd = 15 kHz and I§jhd/Idc = 0-05, while Jf/if, = for other 
cases. 



Figure 4.2 shows several time snapshots of artificial Gaussian eddies (Eq. 
(4.1[>) in the left column and the corresponding synthetic 2D BES data in the right 
column (with DC component removed from Eq. ( 4.5[ >). The eddies are moving 
upward with (vy) = 5.0 km/s. The top left panel in this figure also shows the 
1/e contour lines of the PSFs for the 32 channels (see Figure 3.6 'a)). Snapshots 
for the synthetic 2D BES data are generated with the bandpass frequency filter- 
ing from 10 to 70 kHz to suppress the noise. As the synthetic 2D BES data have 
only 32 spatial points, spatial interpolation is performed using parametric cubic 
convolution technique |84|. 
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Figure 4.2: Left column: four time snapshots of Gaussian eddies, 5n {x, y, t) given 



by Eq. (4.1 1. Right column: the corresponding normalized synthetic 2D BES data 



given by Eq. (4.5) without the DC component. White lines in the top left panel 



show the 1/e contour lines of the PSFs (see Figure 3.6 'a)), and the white asterisks 
show the optical focal points of the 32 channels of the 2D BES system. 
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Part III 

One Further Step in Understanding 
Plasma Turbulence 
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Chapter 5 



Measurements and physical interpretation 



The mean motion of turbulent patterns detected by a two-dimensional (2D) beam 
emission spectroscopy (BES) diagnostic on the Mega Amp Spherical Tokamak 
(MAST) is determined using a cross-correlation time delay (CCTD) method. Sta- 
tistical reliability of the method is studied by means of synthetic data analy- 
sis. The experimental measurements on MAST indicate that the apparent mean 
poloidal motion of the turbulent density patterns in the lab frame arises because 
the longest correlation direction of the patterns (parallel to the local background 
magnetic fields) is not parallel to the direction of the fastest mean plasma flows 
(usually toroidal when strong neutral beam injection is present). This effect is 
particularly pronounced in a spherical tokamak because of the relatively large 
mean rotation and large magnetic pitch angle. The experimental measurements 
are consistent with the mean motion of plasma being toroidal. The sum of all 
other contributions (mean poloidal plasma flow, phase velocity of the density 
patterns in the plasma frame, non-linear effects, etc.) to the apparent mean 
poloidal velocity of the density patterns is found to be negligible. These results 
hold in all investigated L-mode, H-mode and internal transport barrier (ITB) 
discharges. The one exception is a high-poloidal-beta (the ratio of the plasma 
pressure to the poloidal magnetic field energy density) discharge, where a large 
magnetic island exists. In this case BES detects very little motion. This effect is 
currently theoretically unexplained. 




motion 




patterns 



This chapter is largely taken from Ref. Ii66t 



5.1. Introduction 



5.1 Introduction 

It is now widely accepted that turbulent transport in magnetically confined 
fusion plasmas can exceed the irreducible level of neoclassical transport by an or- 
der of magnitude or more dl. However, both theoretical and experimental works 
of the past two decades Il42ti50ll52ll85ti88]l suggest that sheared E x B flows can 
moderate such anomalous transport and hence improve the performance of mag- 
netically confined fusion plasmas. 

With the aim of characterizing the microscale plasma turbulence and search- 
ing for correlations between it and the background plasma characteristics, a 
two-dimensional (8 radial x 4 poloidal channels) beam emission spectroscopy 
(2D BES) system f68^f has been installed on the Mega Amp Spherical Tokamak 
(MAST). It is able to measure density fluctuations at scales above the ion Lar- 
mor radius pi, viz., k±pi < 1, where k± is the wavenumber perpendicular to the 
magnetic field. The 2D BES view plane lies on a radial-poloidal plane at a fixed 
toroidal location. Following the detected turbulent density patterns on this view 
plane allows one to determine their mean velocity in the radial and poloidal di- 
rections. Typically, there are no significant mean plasma flows in the radial direc- 
tion in a tokamak, whereas considerable apparent poloidal motion is detected by 
the 2D BES system. 

In this Chapter, we show experimentally that this apparent poloidal motion 
is primarily due to the strong mean toroidal rotation of the plasma. The fluc- 
tuating density patterns are highly elongated in the parallel direction, and their 
toroidal advection produces apparent poloidal motion due to the projection ef- 
fect. This effect is particularly pronounced in MAST because of a relatively large 
pitch angle of the magnetic field compared to conventional tokamaks. The BES 
measurements are shown to be consistent with a dominantly toroidal mean flow; 
the poloidal flows are of the order of the diamagnetic velocities. These results are 
obtained using the cross-correlation time delay (CCTD) method, which is a fre- 
quently used statistical technique to determine the apparent velocity of density 
patterns |[74ll59ll . We also investigate the method itself thoroughly to determine 
the statistical uncertainties of the technique. This is done by generating synthetic 
2D BES data with random Gaussian density patterns calculated on a graphical 
processing unit (GPU) card using CUDA (Compute Unified Device Architecture) 
programming. 



This Chapter is organized as follows. In Section 5.2 we explain how the ap- 
parent velocity of turbulent density patterns can be inferred from the 2D BES 
data. We also show what physical effects contribute to the apparent velocity cal- 
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culated by the CCTD method. The CCTD method to determine the velocity of 
the density patterns and its statistical reliability are studied in Section 5.3 using 



synthetically generated 2D BES data described in Chapter |4j In Section 5.4 



we 



present the experimental results with the aim of identifying the main cause of 
apparent motion of density patterns measured by the 2D BES system. Our con- 



clusions are presented in section 5.5 



5.2 Velocity of density patterns 

From the time-dependent 2D measurement of density fluctuations, one can 
infer the apparent velocity of the density patterns. This has been the subject of 
much attention Il60l [741 |90l - l94B in the hope that this velocity can be related in a 
more or less straightforward way to the actual plasma flows. We will first ex- 
plore how the mean pattern velocity can be determined and then discuss the 
interpretation of this quantity. 



5.2.1 The cross-correlation time delay (CCTD) method 

The CCTD (cross-correlation time delay) method has been widely used to de- 
termine the apparent velocities of turbulent density patterns detected by BES 
systems, and it is well described in [|74|| and 1189 II . Here, a brief summary of 
the method is provided. The normalized fluctuating intensity of the photons, 
i = 51 /I, measured by a 2D BES system (see Chapter |3| is a function of the 
radial x, vertical (poloidal) y and time t coordinate^ 1 = 1 {x,y,t). The cross- 
correlation function of this fluctuating signal is defined as 

'/ (x, y, t)i{x + Ax, y + Ay,t + At)) 
C {Ax, Ay, At) = ^ ' (5.1) 

P (x, y, t))(p{x + Ax, y + Ay,t + At) 



where Ax and Ay are the radial and vertical (poloidal) channel separation dis- 
tances, respectively. At is the time lag, and (■) denotes time average defined in 
Section 



5.3 



The apparent poloidal velocity of the density patterns detected 
by the 2D BES system can be determined from the time lag At = r^l^y^ at which 
the cross-correlation function reaches its maximum for a given Ay and Ax = Oj^ 
If a straight line is fitted to the experimentally measured r^l^^ {Ay), the inverse 

^The coordinate system used in this Chapter is same as the one used in Chapter |4] 
^We concentrate on the apparent mean 'poloidal' motion of the density patterns. Thus, the 
information about the radial correlations of the 2D BES data is not used in this Chapter. 
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Figure 5.1: (a) Cross-correlation functions calculated using Eq. (5.1)) for Ay = 
0.0 cm (black solid line), 2.0 cm (red dash line), 4.0 cm (blue dash dot line) and 
6.0 cm (green dash dot dot line). t^I^^^ (A?/) is the position of maximum of the 
cross-correlation function, (b) Position of maximum T^l^y. i^y) arid a linear fit. 
The measured velocity is 11.4 ± 0.1 km/ s. 



of its slope is the velocity f^^^. Although any two poloidally separated channels 
are sufficient to determine v^^^, using just two channels is insufficient to estimate 
the uncertainties in the linear fit. Thus, in this Chapter, all four available poloidal 
channels are used to determine these quantities. This assumes that the mean ve- 
locity does not change over the time the density patterns take to move past the 
four poloidal channels and that the lifetime of these patterns is sufficiently long, 
so the same patterns are observed by all four channels. 



Figure 5.1 shows an example of this procedure. This example is based on a 
synthetic data set consisting of Gaussian-shaped random "eddies" moving with 
the poloidal velocity of 10.0 km/ s, which are then used to produce artificial 2D 
BES data (see Chapter |4] for the description of the synthetic data). With the four 
available poloidal channels, cross-correlation functions are calculated using Eq. 



(5.1 1 and shown in Figure 5.1 'a); r^l.^^^ is plotted as a function of Ay in Figure 
5.1 b). The inverse of the slope of a fitted straight line is the velocity v^^'^. Note 



the slight discrepancy between the actual and CCTD-determtned velocities. The 



origin and size of this discrepancy are discussed in Section 5.3 



5.2.2 Physical meaning of the CCTD-determined velocity 

Using the described CCTD method, the 2D BES system on MAST is expected 
to be able to determine Vy^^ as has been done previously on TFTR |i74| and DIII-D 
1 3 using their BES systems BlSHlSl. However, as McKee et al. ||9Zl|98l pointed 
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Figure 5.2: Cartoon illustrating how the mean toroidal plasma flow (Uz) induces 
an apparent mean poloidal motion. An elongated density pattern (shaded oval) 
along the magnetic held line (green dash dot) is advected by the toroidal flow 
(blue arrow). Because the longest correlation direction of the density pattern 
is not in the toroidal direction, the apparent mean poloidal flow (green arrow) 
arises. The apparent velocity is —11^ tan a + Uy —Uz tan a, where a is the local 
magnetic pitch angle. 



out, one must distinguish between the poloidal velocity measured by 2D BES 
system (Vy^^) and the actual velocity of the poloidal plasma flow (Uy). 

The mean plasma flow can be decomposed into toroidal (Uz) and poloidal 
(Uy) components. For typical tokamak plasmas where strong neutral beams are 
injected, \Uz\ ^ \Uy \ is satisfied as any mean poloidal flows are strongly damped 
Il55] - I56ll , leaving Uy of the order of the diamagnetic velocity ~ p*fth/ where p^, = 
Pi/ a, a is the plasma minor radius, and f th is the ion thermal velocity. Note that 
Uz can be on the order of fth for the neutral-beam-heated plasmas. Thus, Uy 
can be ignored compared to Uz, except possibly in regions with strong pressure 
gradients. 

As the 2D BES system on MAST observes the density patterns advected by 
Uz, there will be an apparent motion of the patterns in the poloidal direction, as 



shown in Figure 5.2 This effect is analogous to the apparent up-down motion of 
helical strips of a 'rotating barber-pole' (cf. 162']). The magnitude of this apparent 
velocity can be readily calculated via elementary geometry: namely, we expect 
the BES system to "see", to lowest order in p^, 

vj'^^ ^ -f/^tana, (5.2) 
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where a is the pitch angle of the local magnetic field line. 

Eq. (5.2 1 is experimentally verifiable because all three physical quantities are 
readily obtained by separate diagnostics: f from the 2D BES system, from 
the Charge eXchange Recombination Spectroscopy (CXRS) system [|99lL and a ei- 
ther from EFIT equilibrium reconstruction [|100 | 1 or directly from the Motional 
Stark Effect (MSE) system ilQlllIOl on MAST. Although the CXRS system mea- 
sures the toroidal flow of the C^"*" ions, the difference between the velocity of the 
C^"*" ions and the bulk plasma ions, , is predicted to be on the order of p* in a 
strongly beam-heated plasma | |103 | | . In Section 5.4 Eq. (5.2) will be experimen- 
tally verified for various types of discharges. Agreement will indicate consistency 
of the experiment with the assumptions behind Eq. (5.2). Such agreement will 
indeed be obtained, except in one intriguing case. 

Let us now consider what are the assumptions necessary for Eq. (5.2 1 to hold 
by analysing how the estimated v^^^ depends on actual physical quantities as- 
sociated with plasma flows and fluctuations in a tokamak. The cross-correlation 
function Eq. (5.1 ) of the normalized fluctuating photon intensity / can, in view 
of Eq. ( 3J| ), be considered proportional to the cross-correlation function of the 
relative ion density fluctuation bnjn (by definition, {bn) = 0). Therefore, the 
CCTD-determtned velocity of the density patterns can be related to the actual 
physical quantities in a tokamak by invoking the ion continuity equation. Split- 
ting also the ion velocity into mean and fluctuating parts, u = U + 5u, {6u) = 0, 
we have 

Oti OSti 

— + -— + V ■ (nil + ndu + 6nU + 6n6u) = 0. (5.3) 
ot at 



Averaging this equation and subtracting the averaged equation from Eq. (5.3), 
we obtain 

d5n 

-— = -V ■ (nSu + 6nU + 5n5u - {6n6u)) . (5.4) 
ot 

We will now order various terms in this equation in terms of the small parameter 
P* = Pi /a which is approximately 1/50 — 1/1 00 in MAST. 

Assuming that the spatial scale of all mean quantities is ~ 0{a) while the 
spatial scale of all fluctuating quantities is ~ 0{p^a), and also 5n/n ~ 6u/vth ~ p*, 
we get 

d 5n „5n „ ^ 
= -U-V V-Su 

ot n n 

-5u ■ Vlnn - V ■ (—5u] - — {V-U + U ■ Vlnra) 



n / n 

-O (p^) , (5.5) 
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where we have dropped all terms ~ 0{pl) and smaller. The first two terms on the 
right-hand-side are ~ 0{vth/a) and the following three terms are ~ C(p*fth/a)- 
Note that we have not yet made any assumptions about the nature of the mean 
flow U (beyond it being large-scale) or about time scale of the fluctuations. 

In fact, the p* ordering, which is the standard gyrokinetic ordering | |104[ , can 



take us further. First of all, the mean flow turns out to be purely toroidal to lowest 
order II105I . Operationally, this occurs because of the strong collisional (neoclas- 
sical) damping of mean poloidal flows ||53] - l56| . Thus, U = UzZ + Ui, where z is 
the toroidal direction (locally) and Ui ~ 0(p*) including all poloidal flow^and 
first-order corrections to (radial flows, associated with particle fluxes, are, in 
fact, even smaller). Coupled with the fact that mean quantities have no toroidal 
variation in a tokamak, this means that the fifth term on the right-hand-side of 



Eq. (5.5) is also ~ 0{pl), while the first term can be expressed as 



U-V— = U,- + t/i-V — 

n oz n n 

= -Uz^j^ + -^6- V— + C/i ■ V— , (5.6) 

bzOy n bz n n 

where b = (0, by, bz) is the unit vector in the direction of the magnetic field in a 
local orthogonal Cartesian system {x: radial, y: poloidal and z: toroidal), and we 
have used the identity 6 • V = byd/dy + bzd/dz. Making a further assumption, 
again standard in gyrokinetics, that the parallel spatial scale of the fluctuating 
quantities is ~ 0{a), we conclude that the second term in the second line of Eq. 



(5.6) is Finally, the gyrokinetic ordering also implies that compressibility 



effects (the second term on the right-hand-side of Eq. (53))) are also order p*. This 



is because the lowest-order fluctuating velocity perpendicular to the magnetic 
field is the incompressible E x B drift and the parallel scale of the fluctuations is 
long, so V ■ ~ ^ll'^^ll ~ C> (p*)- 

Combining Eq. ( 53) ) and Eq. (5.6) together with the estimates described 



above, we find 

d 6n d 6n 6n 

^,— + UeffT, =7—, (5.7) 

ot n oy n n 

where U^ff = —Uzby/bz = — f/^tana is the dominant apparent velocity of the 
density patterns {a is the local pitch angle of the magnetic field line). The term 



containing Ueff is the only C(p°) term in Eq. (5.7). The 0(p*) and higher terms 



^Note that the poloidal velocity Uy of the bulk plasma ions has been measured with the CXRS 
system to be only a few km/ s on MAST |106|, which is consistent with Uy ~ ©(p*). Such mea- 
surements are, however, not routinely available for MAST, and one of the goals for this study is 



to confirm that Uy is indeed small. 
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have been assembled in the right-hand-side: by definition, 7 is such that 

7— = --^b-V C/i-V — 

n n n 



This contains, in order of terms, the effects associated with 

(1) parallel variations of the fluctuations, 

(2) mean poloidal flows of bulk plasma ions, 

(3) compressibility of the fluctuations, 

(4) linear response to mean density gradient (drift waves), 

(5) nonlinear effects (turbulence), 

and a slew of higher-order effects of varying degree of obscurity. 



Thus, the right-hand-side of Eq. (5.7) contains all the nontrivial physics of 



waves and turbulence in the plasma. The apparent velocity of the density pat- 
terns detected by the 2D BES system will not be influenced by these effects to 
dominant order — if the orderings assumed above are correct. What it does con- 
tain is the poloidal signature Uej / of the dominant toroidal rotation of the plasma 
— the 'rotating barber-pole' effect discussed at the beginning of this section. In- 
deed, if Eq. ( |5.7[ > holds and its right-hand-side is small, then, to lowest order, 
the density patterns just drift in the y-direction (poloidal) with the velocity Uej /, 



so the maximum of the cross-correlation function Eq. (5.1) will be achieved at 



r = Ay/Ueff. Hence Eq. (5.2 1 for the BES-measured velocity. 

If we are able to confirm Eq. ( |5.2[ ) experimentally, this means that the theo- 
retical considerations employed above are consistent with the experiment. This 
is important because most of the theories of tokamak turbulence rely on such 
considerations. Note that there are no separate diagnostics capable of measuring 



individually all the 0(p*) terms in Eq. (5.8 1. Therefore, the only conclusion one 



can formally draw from Eq. (5.2) holding is that the sum of these terms is small. 



5.3 Assessment of the cross-correlation time delay 
(CCTD) method 



Before we present our experimental results in Section 5.4 let us describe the 

CCTD method and its statistical reliability in more detail. In this section, errors 

involved in determining the mean velocity of the density patterns by the CCTD 

method are examined using the synthetic 2D BES data generated according to 

the procedure explained in Chapter |4j The velocity measured via the correla- 
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tion function (Eq. (5.1 1) is denoted Vy^ and compared with the prescribed value 
(vy) that appears in Eq. (4.4 1, i.e., the mean poloidal velocity of the synthetic 
data. 5.3.1 provides detailed description of the CCTD method used in this work, 
then four types of error are identified for the quantitative comparisons. These 
errors are evaluated in 5.3.3 and 5.3.4 for different values of {vy) and the eddy 
correlation time Tufe. Subsequent sections are devoted to investigating how the 
existence of global (MHD) modes and temporally varying poloidal velocity affect 
the errors. 



5.3.1 Description of the CCTD method 



As defined by Eq. (5.1 1, cross-correlation functions are calculated as time av- 
erages of the data. For a 20 ms-long synthetic data set containing A^totai = 40, 000 
time data points with the sampling time Atgam = 0.5 fis, we want to determine 
^BES ^ time resolution trcs = 1 itT'S- First, a cross-correlation function Eq. (5.1 1 
is calculated on a sub-time window of the synthetic 2D BES data containing Nc 
points, where Nc < tros/Atsam- Then, such cross-correlation functions are aver- 
aged over A^avg consecutive sub-time windows where N^^-^g = (trcs/A^sam)/^c so 
that an averaged cross-correlation function is obtained at every t^es- In this Sec- 
tion, we use Nc = 80, so A'^avg = 25. 

Denoting f{t) and g{t) the time series over a sub-time window from two 
poloidally separated synthetic 2D BES channels, the cross-correlation function 
Eq. (5.1 ) for this sub-time window is: 



Nc 

^sub (^■^^sam) 



Nc-l 



1 

Nc-l 



\ 



Nc-l Nc-l 
k=0 k=0 



(5.9) 



for any integer r with \r\ < Nc — l. Finally, by averaging Csub for Ns,-vg consecutive 
sub-time windows we obtain the smoothed averaged cross-correlation function 
^(rAtsam) from 1 ms-long data points. 

The CCTD method has a serious limitation due to the fact that the sampling 
time Atsam is finite. In order to calculate Vy^^ using only two poloidally sepa- 
rated channels, a line is fitted through two points on a (Ay, r^l^y^) plane as shown 
in Figure 5.1 Id). The first point is located at [Ay, t^I^^^) = (0, 0) by definition, and 
the second point at {Ay,r Atsam)- Then, possible values of Vy^^ are restricted to 
Ay/ (rAtgam) where r is an integer. For the 2D BES system on MAST, using two 

adjacent poloidal channels {Ay = 2.0 cm) with a sampling time Atgam = 0.5 fis, the 
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possible values of Vy^^ are limited to 40.0, 20.0, 13.3, . . . km/s for r = 1, 2, 3, 

Such a limitation may be mitigated by using four poloidally separated channels. 
However, using four channels is not always possible if the channels that are far- 
thest apart are not correlated. To resolve this issue, we use a second-order polyno- 
mial fit on the cross-correlation function C(r Atsam) to locate its global maximum: 
if '^peak is the point where the discrete cross-correlation function C(r Atgam) is max- 
imum, we use the three values of C(rAtsam) at r = rpeak/ '"peak — 1 and rpeak + 1 to 
fit a second-order polynomial. The "true" maximum is found from this fit. We 
denote the time delay at which this maximum is reached by Tp^^^. 



5.3.2 Definition of errors 



For a given set of 20ms-long synthetic 2D BES data, we calculate Vy^^ with the 
time resolution of 1 ms (Section 5.3.1| >. Furthermore, we do this at three different 
radial location^ so that the average of Vy^^, denoted (vy^^), can be calculated 
using 60 values of Vy^^. To make quantitative comparisons between (vy^^) and 



(vy) defined in Eq. (4.4 1, we define four types of error. 
The normalized bias error 



,BES 



Chi: 



i^y) 



(5.10) 



is a quantitative measurement of the systematic discrepancy between the mea- 
sured and the true value. The normalized random error 



errand 



,BES 



,BES 



I (.BES) I 



(5.11) 



quantifies the degree of fluctuation in the measured Vy^^ with respect to (^^^^^)- 



This value may depend on the MHD contribution in Eq. ( 4.5[ > and the temporally 



varying poloidal velocity 6vy{t) in Eq. (4.4). 

Furthermore, as linear fitting is done to determine Vy^^ (see Figure 5.1 1, two 



other types of error are present. The slope of a linear fit can be denoted as Vy ± 
5f fit where Svat is a degree of the uncertainty of the least-square fit|^ Then, the 
normalized mean of Svnt is 

(5.12) 



-fit 



(fafit) 
I (.BES) I' 



*As described in Chapter |4 Vy (t) are identical at all radial locations. One column in the 
middle and t wo co lumns from the edges of the 2D BES channels are used. 

^In Figure 5.1 we plotted t^^^^ as a function of Ay and determined Vy^^ as the inverse of the 
slope of a fitted line. Operationally, we actually plot Ay as a function of t^I^^. so the slope of a 
fitted line is the u™^. 
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and the normalized random error in Svfit is 



a 



.fit 

rand 



I (.BBS) I ■ 



(5.13) 



These two uncertainties together provide an estimation of how well a linear line 
is fitted to given data points. For example, if the assumption that Tufe is long 
enough so that all four poloidally separated channels observe the same eddies 
is not satisfied, then cr^^i^^ becomes large. On the other hand, if this assumption 
is occasionally satisfied, then crf^^^ exhibits such events. Note that error bars of 



the CCTD-determined apparent velocities in Figures 5.11 5.14 in Section 5.4 show 

In the following sections, these four types of error will be evaluated for vari- 
ous values of {vy) and Tnje, and various ranges of Im^hd' fuuD and Svy^^. 



5.3.3 Measuring mean velocity 



To investigate the reliability of the CCTD method described in Section 5.3.1 



for estimating Vy , we generate a number of synthetic 2D BES data sets with 



various values of {vy) while keeping all the other parameters in Eqs. (4.1 1, (4.4 1, 



(4.5 1 and (4.6 1 constant. In real experiments, there is almost always some tempo- 



ral variation of Vy, thus the RMS value of 5vy in Eq. (4.4 1 is set to 5 % of {vy) in 
this subsection. The synthetic 2D BES data are frequency-filtered to suppress the 



noise before the cross-correlation functions are calculated. Figure |53 shows ex- 
amples of (a) Vy (t) generated according to Eq. (4.4 1 with {vy) = 10.0 km/s and (b) 
the original (black) and frequency-filtered (red) autopower spectra of a generated 
synthetic signal. Here, the noise cut-off level is set to be the 5 times the averaged 
autopower level above 900 kHz (green dashed line). 
Figure 



5.4 



-fit 



and (jfjjj^j defined in Section 



5.3.2 



and cal- 



shoWS (Tbias/ 0-rand, CTmean 

culated for values of {vy) ranging from 1 to 100 km/s. The basic conclusions that 
can be made based on these results are as follows: 

(1) For {vy) < 5.0 km/s, the CCTD method is not reliable. This is due to the fact 
that eddies do not live long enough to be detected by all the poloidally separated 
channels. Indeed, it was a priori clear that {vy) < Ay/ru /e could not be measured. 
This translates to {vy) < A.Okm/s for Ay = 6.0cm and r/jje = 15.0 fis, so our results 
are consistent with this simple criterion. 

(2) The CCTD method usually overestimates {vy) (i.e., abias > 0). This can be 
explained by the effective channel separation distance (Ay) being in fact slightly 



less than 2.0 cm because of the overlapping of the PSFs, as shown in Figure 3^ 
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Figure 5.3: (a) Poloidal velocity Vy{t) generated using Eq. (4.4 1 with {vy) = 
10.0 km/ s. (b) Autopower spectra of the original (black) and frequency-filtered 
(red) synthetic BES signals. The green horizontal dashed line shows the noise 
cut-off level, defined to be 5 times the averaged autopower level above 900 kHz, 
and vertical blue dash-dotted lines indicate the low- and high-frequency cutoffs. 
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Figure 5.4: The four types of error defined in Eqs. (5.10l-(5.13l calculated for 
values of {vy) ranging from 1 to 100 km/s. 



(3) The limitation of the CCTD method due to the finite Atgam is successfully 
overcome by fitting a second order polynomial to the cross-correlation function. 



as explained in Section 5.3.1 
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Figure 5.5: Relative discrepancy between the propagation time Tprop = Ay/ {vy) 
and the times t^I^^^ (black) or Tp^^^ (red) at which the cross-correlation function or 
its envelope reaches their peaks for (a) (vy) = 5.0 km/s and (b) 20.0 km/s. 



5.3.4 Effect of the eddy lifetime 



As explained in Section 5.2. 2[ the CCTD method for determining {vy) is based 
on the idea that the peak of the cross-correlation function occurs at t^I^^^ = Tprop, 
where Tprop = ^u/ {vy) is the propagation time of the fluctuating density patterns 
between detectors poloidally separated by the distance Ay. However, r^l^^ will 
not coincide with Tprop if the lifetime tu je of the fluctuations is not long compared 



to Tprop. The failure of the method for {vy) < 5.0 km/s illustrated in Figure 5.4 



is an example of what happens when Tprop is too large. Here, we investigate the 
effect of Tii fe on t^I.^^ quantitatively, via a systematic tu je scan of the synthetic BES 
data. 

Two values (vy) = 5.0 and 20.0 km/s are chosen for this study. For (vy) = 
5.0 km/s, Tprop = 4.0, 8.0 and 12.0 fis with Ay = 2.0, 4.0 and 6.0 cm, respectively; 
for {vy) = 20.0km/ s, they are 1.0, 2.0 and 3.0/is. The peak time r^"^^^ is found using 
the polynomial fitting method described in Section 5.3.1 and (Tprop — ''"peak) /'^prop 
as a function of tu je is plotted for three different values of Ay in Figure 53j It 
shows that r^l^^ underestimates the true Tprop for small values of T/j/e, leading to 
an overestimation of the (vy), consistent with the results shown in Figure 5.4 It is 
encouraging, however, that even relatively low velocities of just a few km/s can 
be determined by the CCTD method with reasonable accuracy (~ 20%). 

It is also possible to consider the global maximum of the envelope of the cross- 
correlation function. We use a Hilbert transform to determine the time delay Tp^^^ 

at which the envelope of the cross-correlation function is maximum [[74| as shown 
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Figure 5.6: Comparison of correlation function and its envelop at Ay = 6.0 cm 
with (vy) = 5.0 km/ s and Tuje = 15.0 ^s. It is clear that the peak position of the 
envelope function r^^l;^ is smaller than that of the correlation function t^I^-^. 



in Figure 5.6 The comparison between Tp^Jj^ and Tprop is shown in Figure 



5.5 



It is 



clear that r^^l]^ has a much stronger dependence on tu je than r'^l^^^, so this measure 
will not be used to estimate {vy) in this work. We note, however, that the strong 
dependence of t^H-^ on the eddies' lifetime Tuff. and of r^l.^^ on their propagation 
time Tprop may provide a way to measure correlation times in the plasma frame. 
Such an investigation is currently being pursued and will be reported elsewhere. 



5.3.5 Effect of coherent MHD modes 

Many experimental 2D BES data sets on MAST exhibit strong MHD (global 
mode) activity in addition to the small-scale turbulence. Removing such global 
modes in the frequency domain is not straightforward as they can have multiple 
harmonics extending into higher frequencies. While they could be filtered out 
relatively easily in the wavenumber domain, constructing wavenumber spectra 
with a very limited number of spatial data points is difficult. Thus, it is useful to 
investigate how the presence of such modes affects the quality of our measure- 
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Figure 5.7: Four types of error (a) as functions of the RMS levels of a global 
mode I Mm) relative to that of turbulence signal 51^^^^; the frequency is fixed 
at fMHD = 15.0 kHz; (b) as functions of the global mode frequency /a/hd at fixed 
Ikihd/^I^^^ = 5-0. Note that atias in (a) is scaled down by a factor of 10, and 
some points are missing in (b) because they are out of the plot range. 



ment of {vy). In this section, this is done by using synthetic BES data sets with 
different RMS levels Im^^d '^^^ frequencies fMHD of the global oscillations (the 
^MHD term in Eq. (4.5)). 



The four errors (abias/ o-rand/ o"mean and a^^^^) are calculated for various ratio of 
^MHD to the RMS value of 6P^ (t) (i.e., 5/^^/5 j^ggg ^j-j-^j-g ^j.^ p^^^. 

ted in Figure 5.7 'a) for the /^hd scan. Here, the frequency of the global mode 
fMHD = 15.0 kHz and (vy) = 10.0 km/s. It is clear that if the power level of the 
mode is larger than that of the turbulence signal, then the CCTD method pro- 
duces large bias errors cfuas- To examine how the frequency of a global mode 
affects the errors, fMHD is varied with a fixed value of Imhd/^^^^^^ = 5.0. The 
results of this scan are shown in Figure 5.7 'b). It shows that (Jbias can be either 
positive or negative with different values of fMHD meaning that global modes 
in real experimental data can cause both over- and under-estimation of the true 

Figure 5.8 shows how different frequencies fMHD can cause such an over- 
or under-estimation of the {vy). Two identical sets of synthetic BES data with 
{vy) = 10.0 km/s are generated, one with and another without a global mode at 
(a) fMHD = 15.0 kHz and (b) fMHD = 40.0 kHz, with = 5.0. Without 

the global modes, the cross-correlation functions with Ay = 6.0 cm (red dashes 
in Figure 5.8 1 have the expected value t^I.^^^ ^ 6.0 fis for both cases. In contrast. 
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Figure 5.8: Cross-correlation functions of the random eddies only (red dash), the 
global mode only (blue dash dot) and the eddies with the global mode (black 
solid) with (a) fuHD = 15.0 and (b) fMHD = 40.0 kHz. Green arrows indicate 
the position of r^l^^^, which does not coincide with the maximum of the cross- 
correlation function of the eddies only (red dash). 



the presence of the global mode in the synthetic BES data shifts r^^^y. towards (a) 
smaller time-lag (over-estimation) or (b) larger time-lag (under-estimation). 

We conclude that a global (MHD) mode with /f/^f, > 61^^^^ affects the struc- 
ture of the cross-correlation functions (both the shape and the position of t^I^^^) 
rendering the CCTD method unreliable. 



5.3.6 Effect of temporally varying poloidal velocity 



No physical quantities are absolutely quiet in real experiments, thus it is nec- 
essary to investigate how the RMS level 5Vy^^^ of the temporal variation of the 
poloidal velocity (see Eq. (4!4)>) influences the measurement of (vy). 



Figure 5.9 shows how finite SVy^^^ / {vy) (with (vy) = 10.0 km/ s) affect the four 



errors defined in Section 5.3.2 It appears that abias saturates at around 50% for the 
scenarios we have investigated, while other three errors increase without show- 
ing any sign of saturation. Thus, the CCTD method to measure {vy) is subject to 
a non-negligible bias error (up to ~ 50%) if the RMS level of temporal variation 
of the poloidal velocity is greater than a half of the mean poloidal velocity. 
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Figure 5.9: Four types of error for various RMS levels (vy) of temporally varying 
poloidal velocities. 



5.4 Experimental results 

In this section, we apply the CCTD method to 2D BES data from MAST dis- 
charges to determine the apparent mean poloidal motion (Vy^^) of the ion den- 
sity patterns. Then, Vy^^ is compared with the 'rotating barber-pole' velocity 
{Uz tan a) where the toroidal plasma velocity Uz is obtained from the CXRS sys- 
tem ||99ll and the local magnetic pitch angle a either from EFIT equilibrium re- 
construction IITOOl or the MSE system IITOTllToa . 

The 2D BES data are first bandpass-filtered from 20.0 to 100.0 kHz to reduce 
the noise level. The low-pass filter removes the high-frequency noise component 
from the photon shot noise and electronic noise as the signal hits the noise level 
above 100.0 kHz in general, while the high-pass filter reduces the contribution to 
the signal from low-frequency, coherent MHD modes as well as the beam noise. 
The apparent mean poloidal velocity of the density patterns Vy^^ is determined 
from average correlation functions calculated over 25 time intervals of 40 /xs dura- 
tion, resulting in total 1 ms averaging. Second-order polynomial fitting is applied 
around the maximum of the correlation function so that r^l^^^ can be obtained 
from a continuous time domain rather than a discretized one due to the 0.5 /is 



sampling time as described in Section 5.3.1 Finally, five consecutive values of 
^BES obtained in this manner are averaged, so the total averaging time is 5 ms 
which is the effective time resolution of Vy^^. Using these five values of Vy^^, the 
time average of various errors defined in Eqs. ( 5.11[ >-( 5l3| > in Section 5.3.2 are also 
computed. 
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We present measurements of Vy^^ from four different discharges: shot #27278 
(L-mode), shot #27276 (H-mode), shot #27269 (ITB) and shot #27385 (high- 
poloidal-beta). All four discharges had double-null diverted (DND) magnetic 
configurations and co-current NBI (neutral beam injection). In all of these dis- 
charges, the 2D BES system viewed at nominal major radial position of i? = 1.2 m 
corresponding to normalized minor radii r/a = 0.2 — 0.3 for L- and H-modes, and 
r/a = 0.3-0.4 for ITB and high-poloidal-beta discharges. The evolution of key pa- 



rameters for these discharges is shown in Figure 5.10 The evolution of plasma 
current, line-integrated electron density and poloidal beta characterize the over- 
all behaviour of plasmas, while the non-zero S-beam voltage corresponds to times 
when the 2D BES system obtains localized density fluctuation. The Da intensity 
trace is used to identify when the H-mode discharge (shot #27276) goes into its 
H-mode: namely, at t = 0.21 - 0.28 s. Note that the ITB discharge (shot #27269) 
starts developing a strong temperature gradient at ~ 0.2 s and the peak ion (C^+ 
from the CXRS) temperature keeps increasing until the NBI cuts off at 0.3 s. The 
viewing position of the 2D BES system is in the middle of the strong temperature 
gradient region for this discharge. 

5.4.1 L-mode (shot #27278), H-mode (shot #27276) and ITB (shot 
#27269) discharges: v^^^ « -Uztana 

Time evolution of (a) cross-power of the fluctuating magnetic field signal from 
two toroidally separated outboard Mirnov coils, (b) cross-power and (c) temporal 
cross-correlation of density fluctuations from two poloidally separated BES chan- 
nels (two mid-channels separated by 2cm) located at -R = 1.21m are shown in Fig- 
ures 5.11 (L-mode), 5.12 (H-mode) and 5.13 (ITB discharge). Here, a cross-power 



is defined as the Fourier transform (in the time domain) of the cross-correlation 
function (Eq. ( |5.1[ >) with finite channel separation. The (minus) apparent mean 
poloidal velocity {—Vy^^, circles) determined by the CCTD method and the 'ro- 
tating barber-pole' velocity {Uz tan a, red solid lines) are also shown in panels (d) 
at i? = 1.13 m and (e) at R = 1.21 m for these three discharges. The error bars rep- 



resent the mean error {6vfit) of the least-squares fit, as discusses in Section 5.3.2 

Despite the fact that these three discharges belong to three very different 
classes, there are common features in the apparent mean poloidal velocity: 
(1) Vy^^ is not reliable (i.e., has large error bars) when strong MHD activity is 
present. The cross-power spectrograms from BES show clear signatures of MHD 
modes with many harmonics, which hamper filtering the BES signal over the fre- 



quency domain. The temporal cross-correlations, i.e.. Figures 5.11 'c), 5.12 'c) and 
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Figure 5.10: Evolution of (a) plasma current (b) line-integrated electron density 
(c) poloidal beta (d) NBI (S-beam) injection energy and (e) edge Da intensity of 
L-mode (shot #27278, black solid), H-mode (shot #27276, red dash), ITB (shot 
#27269, green dash dot) and high-poloidal-beta (shot #27385, blue dash dot dot) 
discharges. 



5.13 'c), also show that these MHD modes have much longer correlation times 



(> 0.3 ms) than the turbulent density patterns. The Mirnov signals show that 
plasmas develop chirping/ fishbone modes with low toroidal and poloidal mode 
numbers [I83II within the filtered frequency band after the shaded region in Fig- 
ures |5^lT|a), |5J2ta) and |5.13[a). The effects of these MHD (global) modes on the 



CCTD method are investigated in Section 5.3.5 where it is found that such activ- 
ity can increase not only the absolute values of the bias errors but also the linear 
fitting errors on Vy^^. Thus, comparisons between Vy^^ and Uz tana are difficult 
to make during the periods where the MHD activity is strong. 
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Figure 5.11: The evolution of shot #27278 (L-mode) showing (a) cross-power spec- 
trogram of the fluctuating B-field signal, (b) cross-power spectrogram and (c) 
cross-correlation of the density fluctuations from BES at i? = 1.21 m. The time 
evolution of —Vy^^ (circles) and the 'rotating barber-pole' velocity {Uz tan a, red 
solid line) at (d) R = 1.13 m and (e) R = 1.21 m. Shaded region shows where 
—yBEs jj^ ^g^^ f^Q^ generate cross-correlation functions and power spec- 
tra are discussed in Appendix |Bj 



(2) During the periods of weak MHD activity, i.e., 0.11-0.15 s for the L- and H- 
mode discharges, and 0.16-0. 22s for the ITB discharge, it is clear that the apparent 
mean poloidal velocity of turbulent density patterns is dominated by the 'rotat- 
ing barber-pole' velocity, i.e., Eq. (5.2) holds, and the sum of all the terms of the 
order of or higher in Eq. (5.8 1 is indeed small. 
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Shot #27276 (H-mode) 



Shot #27269 (ITB discharge) 
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Figure 5.12: Same as Figure 5.11 for 
shot #27276 (H-mode). 



Figure 5.13: Same as Figure 5.11 for 
shot #27269 (ITB discharge). 



Note that the H-mode discharge (shot #27276) goes into its H-mode at ~ 0.21s 
(thus, Vy^^ = —Uz tana is only true before the L-H transition, strictly speaking), 
which can be seen from the intensity trace in Figure 5.10 or from the BES 
cross-power spectrogram in Figure 5.12 the turbulence level drops at the start of 
the H-mode. Any changes of v^^^ during the L-H transition cannot be discussed, 
because the CCTD method with the current data analysis scheme is not reliable 
at this time due to strong MHD activity. 
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Shot #27385 (High-poloidal-beta discharge) 
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Figure 5.14: Same as Figure 5.11 for shot #27385 (high-poloidal-beta discharge). 
Note that neither (b) the cross-power spectrogram nor (c) the temporal cross- 
correlation of the BES signal show any MHD activity. 



5.4.2 High-poloidal-beta discharge (shot #27385): 

Shot #27385 has a relatively higher poloidal beta (the ratio of the plasma pres- 
sure to the poloidal magnetic field energy density) than the three discharges dis- 
cussed in Section [5.4.1 (see Figure 5.10[ >. Thus, it is more susceptible to tearing 
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modes (i.e., formation of magnetic islands) II1071I108 B. The cross-power spectro- 
gram between the two toroidally separated outboard Mimov coils displayed in 
Figure 5.14 'a) shows am/n = 3/2 tearing mode on the g = 1.5 flux surface start- 
ing at ~ 0.11 s; its frequency increases from < 10 kHz to ~ 25 kHz at ~ 0.19 s. 
Then, am/n = 2/1 mode (fundamental frequency < 10 kHz) develops and locks 
to the wall resulting in complete braking of the toroidal rotation of plasmas at 
~ 0.25 s. Here, m and n denote the poloidal and toroidal mode numbers, re- 
spectively, and q for the safety factor. The q = 3/2(2/1) surface moves from 
R = 0.85(0.90) m to 1.24(1.27) m during the initial 0.15 s period and then stays at 
constant R. 

The 2/1 mode is not expected to be seen on the BBS signal as it is bandpass- 
filtered from 20.0 — 100.0 kHz, and no trace of the 3/2 mode is visible in the BBS 
signalj^ Consequently, the Vy^^ determined by the CCTD method does not con- 
tain large error bars during the whole discharge. 

The time evolution of —Vy^^ and Uz tan a in Figure 5.14 'd)-(e) at two different 
radial locations shows that the two velocities do not agree each other at all during 
the period when the 3/2 and 2/1 modes are present. What we find, remarkably, 
is that while the plasma continues to rotate toroidally (as attested by the CXRS 
data), there is virtually no detectable corresponding motion of the density pat- 
terns. In fact, they seem to exhibit a weak rotation in the opposite direction to 
the expected rotating barber-pole effect. Formally, this means that the plasma 
effects in the right-hand-side of Bq. ( 5.7[ > are not small and are able to cancel al- 
most exactly the toroidal rotation, i.e., an effective velocity of the density patterns 
develops in the plasma frame that to lowest order is equal to minus the rotation 
velocity. We do not currently have a theoretical explanation for this effect. There 
is very little apparent difference between the turbulent density patterns in this 
discharge compared to others, except somewhat longer radial correlation lengths. 



5.5 Conclusions 

We have analysed 2D BBS data from different types of discharges on MAST 
to determine the apparent mean poloidal velocities of the ion-scale density pat- 
terns using the cross-correlation time delay method. The dominant cause of the 
apparent poloidal motion of the density patterns is experimentally identified to 
be due to the fact that field aligned patterns are advected by the background, 
dominantly toroidal, plasma rotational flow, i.e., the 'rotating barber-pole' effect 

^Because the mode flattens the mean density profile within the island, shaking of flux surfaces 
does not induce density fluctuations in the BES signal. 
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dominates the apparent mean motion of the density patterns in the lab frame. 
This conclusion holds for the L-, H-mode and ITB discharges we have investi- 
gated. An exception to this rule is found to be the investigated high-poloidal-beta 
discharge, where a large magnetic island is present, and the apparent velocity of 
the density patterns is very small, despite strong toroidal rotation. Identifying 
the causes of this effect by investigating the behaviour of the turbulent density 
patterns quantitatively is left for future work. 
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Chapter 6 



Experimental signatures of critically 
balanced turbulence in MAST 



This chapter is largely taken from Ref. I\109t 

Beam Emission Spectroscopy (BES) measurements of ion-scale density fluctua- 
tions in MAST are used to show that the turbulence correlation time, the drift 
time associated with ion temperature or density gradients, the particle (ion) 
streaming time along the magnetic field and the magnetic drift time are consis- 
tently comparable, suggesting a "critically balanced" turbulence determined by 
the local equilibrium. The resulting scalings of the poloidal and radial correla- 
tion lengths are derived and tested. The nonlinear time inferred from the density 
fluctuations is longer than the other times; its ratio to the correlation time scales 
as i^^^'^^^'^, where z/*, = ion collision rate/ streaming rate. This is consistent with 
turbulent decorrelation being controlled by a zonal component, invisible to the 
BES, with an amplitude exceeding the drift waves' by ~ i^^i^'^- 

6.1 Introduction 

Microscale turbulence hindering energy confinement in magnetically con- 
fined hot plasmas is driven by gradients of equilibrium quantities such as tem- 
perature and density. These gradients give rise to instabilities that inject energy 
into plasma fluctuations ("drift waves") at scales just above the ion Larmor scale. 
The most effective of these is believed to be the ion-temperature-gradient (ITG) 
instability p2l IllOl A turbulent state ensues, giving rise to "anomalous 

transport" of energy [|26|. It is of interest, both for practical considerations of 
improving confinement and for the fundamental understanding of multiscale 



6.1. Introduction 



plasma dynamics, what the structure of this turbulence is and how its ampli- 
tude, scale(s) and resulting transport depend on the equilibrium parameters: ion 
and electron temperatures, density, angular velocity, magnetic geometry, etc. 

Fluctuations in a magnetized toroidal plasma are subject to a number of dis- 
tinct physical effects, which can be thought about in terms of various time scales 
such as the drift times associated with the temperature and density gradients, the 
particle streaming time along the magnetic field as it takes them around the torus 
toroidally and poloidally, the magnetic (V-B and curvature) drift times of parti- 
cles moving across the field, the nonlinear time of the fluctuations being advected 
across the field by the fluctuating E x B velocity, the time between collisions, 
the shear time associated with plasma rotation. Some of these time scales and, 
consequently, the corresponding physics may be irrelevant, while others play a 
crucial role for the saturation of the linearly unstable fluctuations. There has 
been a growing understanding ||112[ , driven largely by theory II113I - I116I , observa- 
tions |117H119| and simulations of magnetohydrodynamic |120i - ll22| and kinetic 
[1141 1123 I I plasma turbulence in space, that if a medium can support parallel (to 
the magnetic field) propagation of waves (and/ or particles) and nonlinear inter- 
actions in the perpendicular direction, the turbulence in such a medium would 
normally be "critically balanced," meaning that the characteristic time scales of 
propagation and nonlinear interaction would be comparable to each other and 
(therefore) to the correlation time of the fluctuations. This means that the turbu- 
lence is not weak and not two-dimensional, unless specially constrained to be so 

cn. 

Beam Emission Spectroscopy (BES) measurements of density fluctuations in 
tokamak plasmas ||68ll69l[73ll96ll98| have made it possible to probe ion-scale tur- 
bulence in these devices directly. In this chapter, we use such measurements 
in MAST, along with the local equilibrium parameters calculated by other diag- 
nostics, to estimate and compare the characteristic time scales of the turbulent 
fluctuations in the energy-containing range. We obtain, for the first time, direct 
evidence that the correlation, drift and parallel streaming time scales are indeed 
comparable across a range of equilibrium parameters (cf. ||1241I125| ') and that the 
magnetic drift time is part of this "grand critical balance" as well. We also find 
indirect evidence that the decorrelation rate of turbulence is controlled by a zonal 
component whose relative importance to the drift-wave-like fluctuations scales 
with the ion collisionality. 

Before presenting this evidence and its implications (e.g., dependence of the 
correlation lengths on equilibrium parameters), let us describe how it was ob- 
tained. 
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6.2 Experimental data and its analysis 

During the 2011 campaign, density fluctuation data from the BES diagnos- 
tic | |68| on MAST were collected in a variety of discharges (including L- and H- 
modes and internal transport barriers). Here we report the data from 39 neutral- 
beam heated "double-null-diverted" discharges, with no pellet injection and no 
resonant magnetic perturbations. The BES system on MAST collects photons 
from a 2D array of 8 radial x 4 vertical locations in the outboard midplane of 
the tokamak, with 2 cm separation between the adjacent channels in either direc- 
tion. The detected photon intensity (mean + fluctuating, / + 51) is used to infer, 
at each location, the density fluctuation level 5n/n = (1//3bes) i^I/I) |69|, where 
/3bes depends on the mean density n and is estimated based on the Hutchinson 
model HTQI (dependence on the mean temperature is weak). As the BES array was 
moved radially for different discharges, our database contains cases with radial 
viewing positions 10 cm < r < 50 cm from the magnetic axis (the minor radius of 
the plasma is ~ 60 cm). 

Local equilibrium parameters are measured by standard diagnostics: mean 
electron densities rie and temperatures by the Thomson scattering system | jl26| , 
impurity ion (C^+) mean temperatures (assumed to equal the bulk ion tempera- 
ture Tj) and toroidal flow velocity by the Charge eXchange Recombination 
Spectroscopy (CXRS) system 1199 L local magnetic pitch angle a by the Motional 
Stark Effect (MSE) system |102|, and further equilibrium magnetic field informa- 
tion is obtained from pressure- and MSE-constrained EFIT equilibria piOOJ. 

We filter the BES data to the frequency interval [20, 100] kHz and calculate the 
spatio-temporal correlation function 

C{Ax,AZ, At) = 

{51 {x, Z, t) 51 {x + Ax,Z + AZ, t + At)) 

where x, Z and t are the radial, vertical and time coordinates, respectively, and 
Ax, AZ and At are the corresponding channel separations and the time lag; (■) 
is the time average over 5 ms periods. At Ax = AZ = 0, the auto-covariances 
{51 (x, Z, t) 51 (x, Z,t + At)) contain not only the physical signal but also photon 
and electronic noise. We remove this effect by applying LED light to the BES 
channels, obtaining 150 different DC levels of BES signal from to 1.5 V, calcu- 
lating the noise auto-covariance (At) at each DC level with the same band 
frequency filter of [20, 100] kHz, then finding Cn (At) whose DC level of the sig- 
nal matches the DC level of the BES data from the MAST discharges, and sub- 

62 



6.2. Experimental data and its analysis 



C(Ax=0, AZ, At) 




-20 -10 10 
&t [|isec] 

Figure 6.1: An example of the cor- 
relation function in the poloidal- 
temporal plane, C (Ax = 0, AZ, At). 
This data was taken at r = 30 cm, 
toroidal rotation speed was f/<^ = 
10 km/ s and magnetic pitch angle 
a = 20°. The direction of maxi- 
mum correlation is the direction of 
the magnetic field (dashed line). 
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Figure 6.2: An example of the corre- 
lation function in the poloidal-radial 
plane, C (Ax, AZ, At = 0). This data 
was taken at the same location and 



the time as in Figure 6.1 



tracting it from the calculated auto-covariances (see Chapter |3] for more detailed 
description). From the correlation function Eq. ( |6.1[ >, we calculate the local char- 
acteristics of the density fluctuations. The spatio-temporal correlation function 



and spatio-spatio correlation function are illustrated in Figure 6.1 and Figure 6.2 
respectively. 

The fluctuation level at each radial location is obtained from the (noise- 
subtracted) auto-covariance function 5n/n = (1//3bes) ^ {5P{x, Z,t)) / 1 at all 32 
locations and then averaged over the four poloidally separated channels at the 
same radial location. 

The correlation length ^y in the direction parallel to the flux surface and per- 
pendicular to the magnetic field is obtained from the vertical (poloidal) corre- 
lation length via ly = Iz cos a, assuming that the parallel correlation length 
is sufficiently long: t\\ ^ ^j^tana. The correlation length tz is estimated us- 
ing four poloidal channels at each radial location (the top channel is the ref- 
erence charmel) by fitting C (Ax = 0, AZ, At = 0) to the function fz {AZ) = 
pz + {1 — Pz) cos [2TiAZ/iz] exp [— \ AZ\ /£z], where pz is a fitting constant that 
serves to account for global structures such as coherent MHD modes (for which 
C (Ax = 0, AZ = oo. At = 0) = Pz 0). In choosing fz (AZ), we assumed wave- 
like fluctuations in the poloidal direction (731 (drift-wave turbulence), with the 
wavelength and correlation length comparable to each other. It is not possi- 
ble to distinguish meaningfully between the two with only four poloidal chan- 
nels. Assuming wave-like structure is essential as in most cases, we find that 

C (Ax = 0, AZ, At = 0) goes negative and/ or is non-monotonic over the vertical 
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extent of the BES array. 

The radial correlation length is estimated using eight radial channels at 
each poloidal location (the fourth channel from the inward side is the reference 
channel). The correlation function C {Ax, AZ = 0, At = 0) is fitted to the function 
fx (Ax) = px + (1 — Px) exp [— I Ax| /ix], where px plays the same role as pz did for 
fz- The values of ix from four poloidal locations are averaged, assuming that the 
radial correlations do not change significantly within the poloidal extent of the 
BES array. Because we have to use the entire array to estimate ix, the number of 
data points for ix is 8 times smaller than for iy. 

To estimate the correlation time r^, we use the fact that the fluctuating density 
patterns are advected poloidally past the BES array with an apparent velocity 
Vy^^ = U^tana due to the toroidal rotation velocity |66 | | (discussed in Chap- 
ter Isl). We fit C [Ax = 0, AZ, At = r^^^^ (^^)) taken at the time delay t^"^^^ (AZ) 
when the correlation function is maximum at a given AZ [|74||, to the function 
fr (AZ) = exp [— I T^l^y^ (AZ) I /tc] . The reliability of this method relies on the tem- 
poral decorrelation dominating over the parallel spatial decorrelation, viz., we 
require Tc <^ i\\ cos a/U^. Anticipating the critical balance assumption Tc ~ i\\/vthi 



||112 | 1 , where Vtu = y'^Ti/rrii is the ion thermal speed, and denoting the Mach 
number Ma = U^/vtu, we estimate that the fractional error in Tc is ~ Ma/ cos a, 
which was never more than 20% in the MAST discharges we used. 

The four quantities 5n/n, iy, ix and Tc are calculated (see Section 3.3 for more 
detailed descriptions) at 8 radial locations (except ix), every 5 ms for all 39 dis- 
charges. All the fits described above are obtained via the mpf it procedure [75J. 
We consider a data point unreliable and remove it from the database if 

(i) / < 0.3 V (the signal-to-noise ratio (SNR) is too low, i.e., SNR< 140 Il68l ). 

(ii) the estimated correlation lengths are smaller than the distance between the 
channels, ix or iy < 2 cm (to be more precise, data points are valid only if ix and 
iy are larger than the size of the point-spread-functions (PSFs) shown in Figure 



3.6 however, most of estimated iy is greater than 10 cm, and 2 cm is a reasonable 



approximation of the PSFs in radial direction), 

(iii) the assumption that plasma rotation is mostly toroidal is suspect, viz., 
l^yBEs _ f/^tana) /v^^^\ > 0.2 (see Ref. ||66l and Chapter |5|), where v^^^ is cal- 
culated at each radial location using the cross-correlation time delay (CCTD) 
method ||74|, 

(iv) the estimated error in the calculation of Vy^^ is > 20% (see Chapter |5] for the 
possible sources of error in the Vy^^), 

{v)pz orpx > 0.5, 

The last two exclusion criteria pick out the cases when MHD modes are too 
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Figure 6.3: (a) Drift time r* = {iy/ Pi) L^/vtu vs. correlation time Tc; (b) 
streaming time r^t = A/tithi = (B/Bp) nr/vthi vs. Tc; (c) magnetic drift 
time tm = (ix/pi) R/vm vs. Tc; (d) perpendicular velocity shear time Tsh = 
[{Bp/B) dU^/dr]]~^ vs. Tc. In all cases, the colour of points represents rji = Ln/ Lt^- 



strong; they are known to degrade the reliability of the BES data II66II . The re- 
maining database contains 448 points. 



6.3 Time scales 



6.3.1 Correlation time vs. drift time 

The turbulence can be driven by radial gradients in the mean ion and electron 
temperatures Tje and density n. Denoting LT^^^ = |VlnTje| and = |Vlnn|, 
the associated time scales are the inverse drift frequencies: 



-1 



Pi Vthi 

i L 



(6.2) 
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where e = '^thi.e/^i.e are the ion (i) and electron (e) Larmor radii, fthi.e = 
■\/2Tj_e/'^j,e the thermal speeds and i^j e = eB/mi^eC the Larmor frequencies. To 
estimate the drift times, we need information about the local equilibrium (Tj e, 
Lt. ^ L„, B) and the correlation length iy, calculated from the poloidal BBS corre- 
lations. 



In Figure 6.3 'a), we compare the drift times with the correlation time Tc calcu- 
lated from the spatio-temporal BBS correlations. We find that = (0.7 ± 0.3)rc, 
where n = min{r*i, r^,„} and the spread is calculated as the root mean square de- 
viation from the mean value. The scaling holds over an order of magnitude in 
either time scale. Thus, the turbulence appears to be driven by the larger of the 
ion temperature or density gradient. However, for Tc < 10 /is, r^e ~ T*i and for 

^ 10 fis, T^n ~ T*i, so we caimot rule out ion-scale electron drive (e.g., trapped 
electron modes [|127 | | or microtearing 11128141301 ). We find no clear correlation of 
T^e with Tc, or with any of the other time scales discussed below. 

6.3.2 Critical balance 

The standard argument behind the critical balance conjecture is causality 



B116B : two distant points on a field line cannot stay correlated if information can- 
not be exchanged between them over a turbulence correlation time. Assuming 
information travels at Vtu, one gets i\\ ~ vtuTc- This cannot be checked directly 
because there are no diagnostics capable of measuring £| on MAST. Consider- 
ing that the inboard side of the torus is a region of "good" (stabilizing) curvature, 
not much turbulence is expected there, so we assume that, at the energy injection 
scale, £| ~ A [I112| , where the distance along the field line that takes a particle 
from the outer to the inner side of the torus is A = nrB / Bp {r is the minor ra- 
dius at the BBS position on the outer side and Bp the poloidal component of the 
magnetic field) |^ Then critical balance means that Tc should be comparable to 



_1 "Lin "im — p -LIU ir '-}\ 

the ion streaming time (the first two equalities are its definition, the last an as- 



sumption). Indeed, we find Tgt = (0.8 ± 0.3)rc (see Figure 6.3 ]b)) 



The balance Tst ~ implies that the poloidal correlation scale is lyj pi ~ A/L*, 



where = min{LT^, L„} |112| . This is tested in Figure 6.4 'a), showing that while 



^As noted above, our method for measuring Tc would instead yield /wthi if Ma > cos a, but 
that would require much stronger rotation (the smallest value in our database is cos a « 0.76). 

^In a conventional tokamak, A « irqR, where q is the safety factor and R major radius, but in 
a spherical tokamak, the local estimate we use is more appropriate. 
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Figure 6.4: (a) Poloidal correlation length iyj Pi vs. kjL^; (b) Radial correlation 



length ^x/Pi vs. kjR. Colour as in Figure 6.3 



the two quantities are certainly of the same order, we do not have enough of a 
range of equilibrium parameters to state conclusively that this theoretically pre- 
dicted scaling works. 

6.3.3 Magnetic drift time and radial correlation scale 

The time scale of the magnetic (V-B and curvature) drifts is 



1 _ Pi^^ ^6 4) 



where we have assumed that the scale length of the background magnetic field is 
R (major radius at the viewing location) and < iy (this will shortly prove cor- 
rect). While magnetic drift physics may matter (in a torus, curvature contributes 
to the ITG drive Il26l ), it does not have to affect scalings, as, for example, it would 
not in a slab and as it did not in the numerical simulations of [|112[ . In contrast. 



Figure 6.3 'c) shows that in the MAST discharges we have analyzed, tm is not 
negligible and scales with t^, similarly to and r^f As tm contains i^, there are 
8 times fewer data points here than in previous two figures, as explained above. 
We find tm = (1-6 ± 0.7)tc. Thus, a "grand critical balance" appears to hold in 
MAST, viz., Tc ~ ~ r^t ~ tm- 

This suggests that the balance of all relevant timescales determines correla- 
tion scales of the turbulence in all three spatial directions. Indeed, balancing 
Tm ~ Tst/ we find the radial correlation scale ix/pi ~ A/i?, the scaling tested 



in Figure 6.4 'b), with a degree of success. This means that the density fluctua- 
tions we are measuring in MAST are not isotropic in the perpendicular plane, 
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Figure 6.5: (a) The nonlinear time associated with density fluctuations, t^^, vs. 
the correlation time t^, (b) their ratio vs. normalized ion collision rate v^i = I'iiTst- 



Colour as in Figure 6.3 



but rather elongated in the poloidal direction iy/ix ~ R/L^ (~ 5 in our data). In- 
terestingly, this clashes with the reported approximate isotropy (ix ~ iy) both in 
Cyclone Base Case simulations p il2| and in measured DIII-D turbulence (where 
^y/^x ~ 1-4 II60II and ix does not appear to depend on Bp II131II ). Whether this is a 
difference between spherical and conventional tokamaks is not as yet clear. 



6.3.4 Nonlinear time 

Since we know the fluctuation amplitude, we can directly estimate the time 
scale associated with the advection of the fluctuations {du± ■ V5n) by the fluctuat- 
ing ExB velocity 5u± = cB x Vy^/i?^. The electrostatic potential Lp is not directly 
measured, but can be estimated assuming Boltzmann response of the electrons: 
5n/n ^ e^p/Tf,. This estimate ignores trapped particles and, more importantly as 
we are about to argue, also does not apply to ion-scale zonal flows (poloidally 
and toroidally symmetric perturbations of with 5n = ||63l |64|). Thus, the 
non-zonal nonlinear time is 

^^Nz^-i _ 1 _ 1 Tj TeCip _ VtMPi TeSn ^^^^ 

i-xiy ^ i-xiy I C Ti Tq i-xiy ^ 



Figure 6.5 'a) shows that r^^^ is always larger than Tc (or the other time scales dis- 
cussed above) and, furthermore, observed to have an inverse rather than direct 
correlation with it. Since turbulence clearly cannot be saturated by linear physics 
alone, this means that our estimate does not capture the correct nonlinear time. 
We conjecture that it is in fact the coupling to the zonal flows, invisible (directly) 
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to BES (because their 5n = 0)|^ that dominates over the nonlinear interaction be- 
tween the drift-wave-like fluctuations represented by r^^ ||34i |48ll63l lT33tiT37 l. 
It has long been suspected that the relative amplitude of the zonal flows com- 
pared to that of the drift waves depends on the ion collisionality [|63l |138i - |140P . 
We can test this expectation by assuming that Tc is the characteristic time associ- 
ated with the coupling of the drift waves to the zonal flows and so depends on 
their amplitude. Figure 6.5 'b) indeed shows a strong collisionality dependence: 
'^c/t^i^ ~ I'^f^^'^, where = uuTst (the ion collision time itself, u^^^, is at least 
an order of magnitude longer than the time scales that participate in the "grand 
critical balance"). We note that a similar scaling is obtained for vs. uuTc and z/j^r^, 
or just straightforwardly for (t^^)~^ vs. uu. t^^ ~ {vthiPi/ix^y)^'^'^^ /Ti, where 
Lp^^ is the amplitude of the zonal potential, this result implies that the ratio of 
zonal to non-zonal component of the turbulence is Lp^^ /ip^"^ ~ i^^i^'^^^'^- A scaling 
popular in theoretical models of zonal-flow-ITG turbulence is z^~^^^ Il63l . 

We note that this situation is qualitatively distinct from what is seen in nu- 
merical simulations of ITG turbulence far from the threshold |112| . where the 
drift-wave nonlinearity appears to dominate (r^^ ~ t^). However, the turbu- 
lence in a real tokamak is likely to be close to marginal and so possibly in the state 
of reduced transport controlled by weakly-coUisionally damped zonal flows ||64| 
and usually associated with the so-called "Dimits upshift" of the stiff-transport 
threshold ll34l [135l[14n[T42l . 



6.4 Discussion and conclusions 

Our results support the notion that the statistics of turbulence are determined 
by the local equilibrium properties of the plasma. We find little correlation be- 
tween the quantities reported above and the radial locatiori|^ (note that we have 
limited our consideration to temporal and spatial scales and did not touch on 
the fluctuation amplitudes or transport properties, which do of course depend 
on radius). Our results also appeared insensitive to (i.e., not measurably corre- 
lated with) three other parameters that might in principle have proven important: 
Tj/Tg (varied between 0.5 and 2), the magnetic shear s = dlnq/dlnr (varied be- 
tween —1 and 5) and the perpendicular component of the toroidal velocity shear 



Tgj^ = {Bp/B)dU^/dr. In much of our data, r^h > r^Tst (see Figure 6.3 'd)), so it 



■'Zonal flows can be detected using the BES system by looking at the velocity modulation of 
5n II132I which requires more sophisticated statistical analysis. However, this does not change 



our estimation of r^^^. 



^There is a slight bias in Figure 6.5 for larger i^^i to be found farther from the magnetic axis. 
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stands to reason that the statistics of the turbulence would not be dramatically 
affected; in the instances of Tsh ~ Tst, the effect of r^h could not be isolated. In 
general, we expect that a strong velocity shear would change £| via a modified 
critical balance: if r^h < Tst, then ~ ~ i\\ /v^u, so l\\ ~ f thjTsh < A. It would be 
interesting to investigate higher-rotation plasmas, as t~^, when sufficiently large, 
is expected to have a dramatic effect on transport ||4T] - l44ll48ll50ll88lll43| : even in 
our database, there is in fact some evidence that velocity shear might raise the 
critical temperature gradients [|144[ , but we see no signature of this effect in the 
correlation properties of the turbulence. 

We have presented experimental results statistically consistent with a turbu- 
lent state in MAST set by the local equilibrium and in which the time scales of the 
linear drive, turbulence decorrelation, ion streaming and magnetic drifts are all 
similar and scale together as equilibrium parameters are varied. This "grand crit- 
ical balance" implies a three-dimensionally anisotropic turbulence, with parallel, 
poloidal and radial correlation lengths having different parameter dependences 
and £|| ^ ly > ix- Our results also suggest the presence of a zonal component 
with an amplitude i^^^'^^^'^ greater than the drift-wave density fluctuations. Fur- 
thermore, it provides a way to estimate the turbulent fluxes (at least an order of 
magnitude) solely based on equilibrium parameters as we can predict the spatial 
and temporal structures of turbulence and its fluctuation level. 
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Better Tokamak Performance 
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Chapter 7 



Local dependence of ion temperature 
gradient on magnetic configuration, 
rotational shear and turbulent heat flux in 

MAST 



This chapter is largely taken from Ref. I\144t 

Experimental data from the Mega Amp Spherical Tokamak (MAST) is used to 
show that the inverse gradient scale length of the ion temperature R/Lt^ (normal- 
ized to the major radius R) has its strongest local correlation with the rotational 
shear and the pitch angle of the magnetic field (or, equivalently, an inverse corre- 
lation with q/e, the safety factor/the inverse aspect ratio). Furthermore, R/Lt^ is 
found to be inversely correlated with the gyro-Bohm-normalized local turbulent 
heat flux estimated from the density fluctuation level measured using a 2D Beam 
Emission Spectroscopy (BES) diagnostic. These results can be explained in terms 
of the conjecture that the turbulent system adjusts to keep R/LTi close to a cer- 
tain critical value (marginal for the excitation of turbulence) determined by local 
equilibrium parameters (although not necessarily by linear stability). 

7.1 Introduction 

A key physics challenge posed by magnetically confined plasmas in fusion de- 
vices is how the internal energy can be kept from being transported too fast from 
the core to the periphery. This problem is primarily one of turbulent transport, 
the temperature gradient between the edge and the core of a toroidal plasma 



7.1. Introduction 



supplying the source of free energy for the excitation of the turbulent fluctua- 
tions that then act to enhance the effective thermal conductivity and relax the 
gradient. It is expected that it is the ion temperature gradient (ITG) that gives 
rise to the most virulent instabilities (on ion Larmor scales) Il32llll0lllll[ and that 
this is then self -consistently limited by the resulting turbulence ||26|. If we view 
the edge ion temperature as fixed by the physics and engineering aspects of the 
tokamak design that will not concern us here II145L the key question is how to 
maximize the ion temperature gradient. We therefore wish to inquire, experi- 
mentally, on what this gradient depends and how. Motivated by the fact (or the 
conjecture) that the state of the ion-scale microturbulence is largely determined 
by the local (to a given flux surface) equilibrium conditions 111051 11091 1146I - I148I 
and in turn acts back to adjust them locally, we ask what local parameters are 
most strongly correlated with the corresponding value of R/Lt^, the inverse ra- 
dial gradient scale length of the ion temperature {L^^ = \d In Ti/dr\) normalized 
to the major radius R of the torus. 

One may wonder how universal any such measured dependences are likely 
to be for situations with different global conditions, e.g., different neutral-beam- 
injection (NBI) heating powers. It has been recognized for some time that the 
turbulent heat flux tends to increase very strongly (much faster than linearly) 
with R/Lt^ — a phenomenon known as "stiff" transport |l33l|49H50HI12l[ll91. If 
(or when) the transport is indeed stiff, any experimentally measured relationship 
between R/Lxi and other equilibrium parameters should be quite close to some 
critical manifold in the parameter space separating dominant turbulent transport 
from an essentially non-turbulent or weakly turbulent state ("the zero-turbulence 
manifold" |411). This critical manifold would be independent of the power input 
and can be thought of as a local parameter dependence of the critical tempera- 
ture gradient R/Lt^,c = /(g, s, U'^, R/L^, R/Lt,, va, Ti/T^, /Si,... ), where q is the 
safety factor (number of toroidal revolutions per one poloidal revolution of the 
magnetic field around the torus on a given flux surface), e = r/R the inverse as- 
pect ratio (r is the minor radius of the flux surface), s = d\nq/d\nr the magnetic 
shear, U'^ = dU^p/dr the radial shear of the mean toroidal rotation velocity t/^, L„ 
and Lt^ the gradient scale lengths of the plasma density and electron tempera- 
ture, the ion collision frequency, Tj/Tg the ion-to-electron temperature ratio, 
/3j = STrnTi/B"^ the ion-to-magnetic pressure ratio and ". . . " stand for everything 
else (e.g., the many parameters required to fully describe the magnetic configu- 
ration)!^ An important caveat is that the "zero-turbulence" threshold need not be 

^We do not suggest causality between all these parameters and R/Lxi — all local equilibrium 
characteristics, including R/Lt^, jointly adjust to form the critical manifold. 

73 



7.1. Introduction 



the same as the threshold for the existence of linearly unstable eigenmodes. Two 
known examples when this is not the case are the so-called "Dimits upshift" of 
above the linear stability threshold Il34lll35iri42ll and the case of sufficiently 
large U'^ when the system may be linearly stable but strong transient excitations 
p9ll40| lead to sustained subcritical turbulence ||42ll43ll85| . Thus, in general, there 
is a nonlinear threshold with some definite dependence on local equilibrium pa- 
rameters. 

Recent theoretical ||39ll40| and numerical ||4T| investigations suggest that q/e 
and U'^ may be the most important such parameters, at least at low s. Let us 
explain why this is. It is well known, both from experimental measurements 
Il66l [10611 and theory |l55l |56l [lOl, that strong (finite-Mach) flows in a tokamak 
are predominantly toroidal (certainly when plasma is heated by tangential neu- 
tral beams, which produce a toroidal torque)]^ Therefore, unless the magnetic 
field is purely toroidal, any radial shear in the toroidal flow results in sheared 
flow in both the perpendicular {U'j_ = {Bg/ B) U'^, Be is the poloidal field) and 
parallel (f/|| = {B^/B)U'^, B^ is the toroidal field) directions. While perpen- 
dicular flow shear is known (theoretically ||42] - |48l 11501 1151II and experimentally 
| |49]452l 11491 1152[ ) to suppress turbulence and the associated transport, parallel 
flow shear can drive turbulence via the "parallel-velocity-gradient" (PVG) insta- 
bility ||38] - |40B . The average ratio of these two shearing rates on a flux surface, 
f/||/f/^ = B^/Be, can be approximated by q/e and so the degree to which sheared 
equilibrium flow suppresses or drives turbulence is expected to depend on this 
parameter. Indeed, numerical studies of ITG- and PVG-driven turbulence have 
shown that the critical R/Lt^^c at any given value of U'^ increases with decreasing 
q/e {at least for low s t4l|l): while at any given q/e, R/Lt^^c increases with increas- 
ing U'^ provided the latter is not too large [|4T] - |43ll85l (as in most real tokamaks). 

A comprehensive numerical parameter scan of the dependence of R/Lt^^c on 
all other potentially important local quantities (s, R/Ln, R/Lt^, ua, Ti/Te, etc.) 
is probably unaffordable in the near future, so faster progress can be made experi- 
mentally. In this chapter, our first goal is to establish, based on a relatively sizable 
dataset for MAST, what the most important parameters for the critical manifold 
are: we will show that, indeed, the local value of R/Lt^ is most strongly corre- 
lated inversely with the local q/e and positively with the local rotational shear — 
consistently with the result obtained in [41]. 

Our second goal is to obtain an experimental signature that the measured 
R/Lt^ is determined by — or, more precisely, correlated with — the local charac- 

^Any mean poloidal flow exceeding the diamagnetic velocity would be damped by collisions. 
15311511. 
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teristics of the ion-scale turbulence, directly measured by the 2D beam emission 
spectroscopy (BES) diagnostic ||68|. We will show that not only does a strong cor- 
relation between R/ Lt^ and an estimated turbulent heat flux level exist but its (at 
the first glance, counterintuitive) inverse nature is consistent with R/Lx^ staying 
close to the critical threshold R/ Lt^^c and hence with stiff transport. 



7.2 Equilibrium parameters 

A database was compiled of equilibrium quantities (and turbulence character- 
istics; see below) from 39 neutral-beam-heated discharges from the 2011 MAST 
experimental campaign. These discharges had a double-null diverted (DND) 
magnetic configuration, no pellet injection and no applied resonant magnetic 
perturbations. Mean electron density rie and temperature Te were measured with 
the Thomson scattering system | |126 | | , mean impurity ion (C^+) temperature Tj 
and the toroidal flow velocity with the Charge eXchange Recombination Spec- 
troscopy (CXRS) system [99J (we assumed that in these discharges the impurity 
and bulk ions have negligible differences in their temperature and flow veloci- 
ties |153| V The local magnetic pitch angle (Bg/B^) was measured with the Mo- 



tional Stark Effect (MSE) diagnostic II102II ; pressure- and MSE-constrained EFIT 
equilibria [|100| were used to obtain the field strength B. All parameters were 
determined over 5 ms intervals either by averaging if the diagnostic's temporal 
resolution was smaller or by interpolation if it was larger than 5 ms. Only data 
points from a limited range of minor radii 0.6 < r/a < 0.7 (r = a is the edge 
of the plasma) were used, in order to minimize any correlations between vari- 
ous quantities due to their profile dependence alone (thus, we did not attempt to 
prove locality here; see, however, |109I). In total, 988 data points were available. 

From this information, we constructed 7 local dimensionless parameters, 
which, motivated by theoretical models or common sense, we deemed a priori 
the most important ones (we also give the range of variation of each parameter): 
R/Lt, e [0.08,20.3], q/s E [4.0,16.3], s E [1.2,6.0], 7e = U'^t,, = -KrU'^v^u e 
[0.005, 2.5], E [0.04, 13.8], /?/Lt, E [1.43,22.7], u^i = u^^T,t E [0.003,0.12], 
Tj/Te E [0.5, 1.7]. The ion collision rate uu and the perpendicular velocity shear 
U'j_ (which is used instead of U'^) were normalized, as in B109I , to the ion parallel 
streaming time Tgt = A/vthi, where Vthi = ^/2Ti/mi and A = nrB/Bg is the con- 
nection length (the approximate distance along the field line from the outboard to 
the inboard side of the torus, expected to determine the parallel correlation scale 



of the turbulence |112| : if the flux surfaces had been circular, A ^ nqR). The local 
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magnetic configuration is represented by only two parameters: q/e and s. The 
choice of q/e was motivated by the physical considerations outlined in Section 



7.1 (since e varied very little in our database, we do not claim to distinguish any 
individual correlations of R/Lt^ with q and e). It is left for further study whether 
other properties of the flux surfaces matter (e.g., Shafranov shift, triangularity, 
elongation, etc.; some of these may, in fact, affect the stiff-transport threshold 
[1421 1154II ). We have not included n, Tj, Te, which are not normalizable by any 
natural local quantities; note that R/Lt^ usually has a large but trivial correlation 
with Ti. larger temperature gradients lead to larger temperatures in the core. We 
also have not included (3i = SvmTj/i?^ because, in the absence of large variation 
of i? in our dataset, /3j is simply the normalized ion pressure and, similarly to Tj, 
has a large positive correlation with R/Lt^ (it remains to be investigated whether 
the larger level of magnetic fluctuations at larger (3i is large enough to have a 
nontrivial effect on turbulent transport lfT29l[T30l[T55tiT59l ). 

7.2.1 Correlation analysis 

We perform a Canonical Correlation Analysis (CCA) |160| with R/Lt^ treated 



as the dependent variable and the other 7 local parameters itemized above as 
independent ones. This amounts to finding the maximum correlations between 
\iii{R/LtJ and linear combinations of logarithms of 1, 2, 3, . . . , or 7 other parame- 
ters, leading to an effective statistical dependence 

" ^'^^■s?cf#r*.«f#v'm"'. (7.1) 



This is of course not valid if the dependence of R/Lt^ on any of the parameters 
is non-monotonic. A non-monotonic dependence on 7^; is, in fact, expected, with 
R/Lxi first increasing, then decreasing at larger values of due to increased 
transport from the PVG-driven turbulence [|4l] - l43l |85| . However, the range of 
values of 'Je in our database does not extend to sufficiently high values for such 



a dependence to be observed (see Figure 7.1 1. 



The results are shown in Table |7.1t where the values of the canonical corre- 
lation (i.e., the correlation coefficient between the logarithms of R/Lt^ and the 
right-hand of Eq. ( |7.1[ >) are given together with the corresponding exponents 
07. We start by calculating the individual correlations of R/Lt^ with each 
of the 7 parameters and then include pairs, triplets, etc., only if the correlation 
improves. We see that the strongest individual correlation of -R/^Ti are with q/e 
(61%) and 'Je (46%). The overall fit is measurably improved (66%) if both are 

included. Including further parameters does not make a significant difference; 
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Table 7.1: Results of CCA performed assuming Eq. (7.1). Wherever appears, 
that means that the CCA was performed without including the corresponding 
parameter. 
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the third strongest (although not very strong) dependence is on z/*j. The depen- 



dence of R/Lxi on q/e and 7^; is shown in Figure 7.1 which generally increases 
with decreasing q/e and increasing Je^I This is broadly consistent with the ex- 



pectations based on intuitive physical reasoning (explained in Section 7.1) and 
on the numerical study of |j4T|. We note that s was set to zero in the numerical 
study reported in f^] which is different from our experimental cases. However, 
the values of s will not change the quantitative trend of observed R/ Lt^ as a func- 
tion of q/e and as reported in ||42| with s = 0.8. 

The conclusion is that, at least on a very rough qualitative level, it is sensible 
to consider R/Lt, to be a function primarily of q/e and ^e- We note that equi- 



librium database shown in Figure contains time periods when MHD modes 



are active unlike the BES database (see Figure 7.3 1. As a result, it is possible that 
our database may include some outliers due to a fast change of heat transport 
during the MHD activities. Since the q profile tends to change more slowly in 
tokamaks than other equilibrium profiles]^ it may be useful to think of a critical 
curve R/ LTi,c{lE) P1431 parametrized by q/e |j4T|, the latter quantity containing 
the essential information about the nature of the magnetic cage confining the 
plasma. 



^The broad scatter of data points in Figure 7.1 a) suggests that the correlation between q/e 



and 7b is weak; the lack of higher values of 7^; at large q/e is due to the fact that the flow shear is 
weak at earlier times in the discharges, when the central value of q is high. 

*It can be proven that the functional dependence q{ip), where il) is the flux-surface label, only 
changes on the resistive timescale of the mean magnetic field II159I . 
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Figure 7.1: The dependence of RjLx^ (color) on qje and = T^rU'^/vthi, showing 
(a) the individual data points and (b) mean values of R/Lxi within rectangular 
bins. 



7.2.2 CoUisionality dependence 

Even though the z/*j dependence of R/Lt^ is not as strong as q/e and ^e, some 
discernible inverse correlation between R/ and v^i might be expected because 
zonal flows, believed to suppress turbulence Il63l 11351 , should be more strongly 
damped at higher ion coUisionality | |138I4141 |. To isolate this dependence, we 
selected data points for approximately fixed 7^ e [0.7,0.8] and q/e E [5,6] (the 
largest number of data points could be found within these narrow ranges and no 
measurable correlation between R/Lxi and q/e or ^e was present). The resulting 
Figure 7.2 confirms a degree of inverse correlation between R/ Lt^ and v^i. 



7.3 Turbulent heat flux 

The turbulent ion heat flux through a given flux surface is (very approxi- 
mately!) Qi ~ nTiXi/Lxi, where the effective turbulent diffusivity is Xt ~ Su^Tc, 

6u ~ {c/B)ip/iy is the (radial) fluctuating E x B velocity, Tc its correlation time, iy 
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Figure 7.2: R/Lt^ vs. v^i for fixed e [0.7, 0.8] and q/e E [5, 6]. 



its poloidal correlation scale and (f the fluctuating electrostatic potential. The lat- 
ter can be estimated from density fluctuations using the approximation of Boltz- 
mann electrons: e^p/Te ~ 5n/n (e is the proton charge, n and 5n the mean and 
fluctuating density, respectively). Both theory of ITG turbulence | |112 | | and the 
BBS measurements in MAST II109I suggest that ~ r^i = lyL^JviuPi (the drift 
time; pi is the ion Larmor radius). Collecting all this together, we estimate the 
gyro-Bohm-normalized turbulent ion heat flux: 

Qi,turh Pi f RV fTeSnY _ - 

QgB 4 \PiJ \Ti n J 

where Q^b = nTiVtupf/R^- 

Since ion-scale density fluctuations in MAST can be measured directly by the 
EES system, Qturb can be obtained independently of any transport reconstruction 
models such as TRANSP B161L The method of determining 5n/n and ^y using 
the BES system on MAST (8 radial x 4 vertical channels with spatial resolution of 
~ 2 cm II68II ) is explained in detail in II109I . This is done from the covariance and 
correlation functions of the photon intensity fluctuations, averaged over the same 
5 ms intervals for the same 39 discharges as the equilibrium quantities studied 
above, although not in all intervals there was good BES data. Restricted to the 
radial range 0.6 < r/a < 0.7, the number of available data points for 5n/n and iy 
was 102. 



7,3.1 Inverse correlation between R/Lt^ and Q 



turb 



It is shown in Figure 7.3 'a) (Qturb vs. q/e and ^e'> cf. Figure 7.1 1 and Figure 
7.3'h) (Qturb vs. R/Lt/) that smaller Qturb is observed where R/Lt^ is large and 



vice versa. This is perhaps counterintuitive as one might expect that the larger 
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Figure 7.3: (a) Qtu rb (c olor) calculated from BBS dat a ac cording to Eq. (7.2) vs. q/e 
and 'jE (cf. Figure [TJ]); (b) Qturb vs. R/Lt,{'^L Figure [7!4|; (c) Qturb (color) vs. R/Lt, 
and for a fixed range oiq/e E [9, 11], indicated by the dotted horizontal lines in 
(a); open squares are data from Figure [7J 'a) for which BES measurements were 
not available (cf. Figure 7.5 1. 



R/Lt^, the more turbulent the plasma and so the larger the turbulent heat flux. 
This would indeed have been the case had R/Lt^ been externally fixed, as in lo- 
cal flux-tube simulations, where Qturb does increase with R/Lt^ p4llll2| . In con- 
trast, in the real plasma, a certain amount of power flows through a flux surface 
and, if transport is stiff, the temperature gradient (along with other equilibrium 
quantities) adjust accordingly, to stay close to the critical gradient defined by the 
manifold R/ LT^^dlE, q/^)- Indeed, in plasmas with both power and momentum 
injection, there is a regime with R/Lt^ close to R/Lt^ c where the turbulent and 
the neo-classical (collisional) transport are comparable and in which a larger heat 
flux results in lower R/ ■ For a detailed explanation, we refer the reader to |143| 
(see also Fig. 4(b) of [43J (or Figure 7.4[ >, which should be compared to our Fig- 



ure 



7.3 'b)). In brief, in the neoclassical regime, the momentum transport is much 



less efficient than the heat transport, while in the turbulent regime, they are com- 
parable (the turbulent Prandtl number is order unity [|42l |85l 11511 11621 while the 
neoclassical one is small Il55l ), so, as a larger heat flux takes the system (slightly) 
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II/Q = (Data) 
(b) n/Q = 0.04 (Int.) 

n/Q = 0.06 (Int.) 
n/Q = 0.08 (Int.) 
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Q = 2.6 
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Figure 7.4: Numerically generated 
total heat flux Q vs. R/Lt^ show- 
ing that there exist regions where 
a smaller Q corresponds to a larger 
R/Lxi- Figure is taken from Fig. 4(b) 
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Figure 7.5: Analytically calculated 
K = R/Lt^ vs. -fE with Qi > (^2 > 
Qa. Point A3 has a larger R/Lt^ 
than ^2 despite < Q2- The red 
line is the critical temperature gra- 
dient R/ L^i^c below which there is 
no turbulence. Figure is taken from 



Fig. 3(b) of IIT431I 



farther from the marginal state, this leads to much more efficient momentum 
transport, hence smaller velocity shear 7^;, hence a regime with less suppression 
of turbulence and smaller R/Lt^^c (see Figure 7.3 'c), where the correspondence 
between larger Qturb, lower R/Lt^ and lower 7^ is shown at approximately fixed 
q/e G [9, 11]; cf. Fig. 3(b) of IIT43B (or Figure |7!5|)). We stress that all of this happens 
quite close to marginality and so our experimental observation of an inverse cor- 
relation between R/Lxi and Qturb provides circumstantial evidence that R/Lt^ in 
MAST is indeed close to its critical valuejf] 



7.4 Conclusion 

We have found that the normalized inverse ion-temperature-gradient scale 
length R/ has its strongest local correlation with q/e and the shear in the 
equilibrium toroidal flow: R/Lt^ increases with increasing shear, which is a well 
known effect, and with decreasing q/e, which corresponds to an increasing ra- 
tio of the perpendicular to the parallel shearing rates. We note that a similar 
dependence of R/Lt^ (g/e, 7^;) is also observed in JET B163P , suggesting that the 
inverse correlation between R/Lt^ and q/e is perhaps ubiquitous, as would be 
the case if R/ Lt^ were generally fixed at some locally determined critical value 
||4T | . Furthermore, we have found an inverse correlation between R/Lt^ and the 

^This situation is distinct from the experiments on transport stiffness on JET [|49l l50l [1491 ir* 
that they provided vigorous extra heating power using localized ion-cyclotron-resonance heating 
(ICRH) to depart far from marginality. 
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gyro-Bohm-normalized turbulent heat flux (estimated via direct measurements 
of density fluctuations) and argued that this is consistent with R/Lt^ always re- 
maining close to a critical manifold R/Lt^^c {q/^i 1e) separating the turbulent and 
non-turbulent regimes 1431 185111431 (stiff transport). It is thus plausible that we 
have essentially produced this critical manifold for the MAST discharges we in- 
vestigated. Practically, our results suggest that R/Lt^ can be increased by low- 
ering the q/e, which is relatively easier and less expensive than increasing the 
shearing rate in tokamak operations n 



^With tangential NBI heating, it is difficult to increase the toroidal Mach number — and hence 
the equilibrium flow shear — because of the fixed ratio of injected torque to power at a fixed 
injection energy. 
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Chapter 8 
Conclusions 



The poet should say what sorts of thing might happen, that is, the things 

possible according to likelihood or necessity. 

- Aristotles, Poetics Ch. 9 

In this work, we have presented experimental results on the characteristics of 
plasma turbulence measured by the 2D beam emission spectroscopy (BES) sys- 
tem on the MAST spherical tokamak. 

First, we have described the principle of the 2D BES system used to mea- 
sure ion-Larmor scale fluctuating densities up to a few 100 kHz down to a few 
0.1 % level. We find from the generated point spread functions (PSFs) that the 
radial spatial smearing does not depend significantly on the observation loca- 
tions; whereas the poloidal spatial smearing becomes greater as the observation 
locations are moved toward the edge of MAST due to increasing magnetic field 
pitch angles. The PSFs at different major radii vary which complicates the decon- 
volution of the measured signal. However, as the measured poloidal correlation 
lengths are found to be larger than the poloidal widths of the PSFs, the deconvo- 
lution of the signal would only slightly correct the measured correlation lengths. 
BES data with shorter radial correlation lengths than the radial widths of the 
PSFs are not used so that our results are not biased even if the deconvolution is 
not performed. Then, detailed procedures of obtaining statistical properties of 
turbulence such as spatial and temporal correlation lengths and generating syn- 
thetic 2D BES data are discussed. 

The first physics result we have obtained from the measurements is on the 
subject of interpretation of mean motion of fluctuating density patterns. Using 
the cross-correlation time delay method whose statistical reliability is extensively 
investigated with the synthetic 2D BES data, the apparent poloidal motion of 
fluctuating density patterns in the lab frame is found to arise because the longest 



correlation lengths of the patterns are not parallel to the dominant toroidal flow 
of the bulk plasma. This projection effect holds for the investigated L-, H-mode 
and internal transport barrier discharges. An interesting exception to this rule is 
found for the investigated high-poloidal-beta discharge, where a large magnetic 
island exists. Whether this exception is related to presence of the island is not 
known, nor do we successfully identify any feasible physical explanations for 
this phenomenon. Further investigation is left for future work. 

With the understanding of mean motion of turbulence, spatial and temporal 
characteristics of it are compared with the local equilibrium quantities. Such com- 
parisons are performed in terms of various time scales each of which is related 
to a distinct physical effect. The measured turbulence is shown to be consistent 
with the idea of turbulence being critically balanced, that is, the correlation time 
of the turbulence, the drift time associated with the background ion-temperature 
or density gradients, the ion streaming time along the magnetic field line and the 
magnetic drift time are consistently comparable. The balance between the drift 
time associated with the background gradient and the ion parallel streaming time 
provides how the perpendicular (with respect to the magnetic field within a flux 
surface) correlation length of turbulence scales with the ratio of the connection 
length to the gradient scale length; while the balance between the magnetic drift 
time and the ion parallel streaming time gives the scaling of the radial correlation 
length as a function of the ratio of the connection length to the size of a tokamak. 
Furthermore, estimation of non-zonal component of nonlinear time of the turbu- 
lence seems to indicate that the measured turbulence is dominantly decorrelated 
by the scattering of the drift waves by the zonal flows rather than the nonlin- 
ear interactions between the drift-wave-like fluctuations themselves. It is found 
that the ratio of zonal to non-zonal components of the turbulence scales inversely 
with the ion collisionality. 

Finally, we have experimentally shown that the normalised inverse ion- 
temperature-gradient scale length R/Lt^ has the strongest local correlation with 
q/e (safety factor /inverse aspect ratio) and the shear in the equilibrium toroidal 
flow. Furthermore, we argued that the counterintuitive observation of the in- 
verse correlation between the gyro-Bohm-normalised turbulent ion heat flux and 
the R/Lt^ is evidence that observed turbulence level is quite close to its critical 
value indicating that we have essentially produced the critical manifold in the lo- 
cal equilibrium parameter space for the MAST discharges. In practice, we can use 
this extra experimental knob of g/e to control R/Lt^ which is much cheaper and 
easier than the 'shear flow' knob. Thus, the q/e may well become a key parameter 
to ignite the plasmas in the ITER. 
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Appendix A 
Drift waves 



In this chapter, we consider three types of drift waves associated with gradi- 
ents of density, ion temperature and parallel velocity. First, we provide a pictorial 
description of a density-gradient driven drift wave for easier understanding of 
physical mechanism of the wave following ref. | j26| . Then, we use fluid equa- 
tions, being more quantitative, to describe ion-temperature-gradient (ITG) and 
parallel-velocity-gradient (PVG) driven drift waves following Cowley's lecturej^ 

In all cases, we assume that perturbations are electrostatic for the simplicity, 
and we use a cartesian coordinate system where x is the inhomogeneous direc- 
tion, i.e., equilibrium quantities such as density, temperatures and plasma flows 
are varying in the x— direction, the y-direction is taken to be periodic, and the 
z-direction is parallel to the homogeneous straight magnetic field lines, B = Bz. 
Note that the assumption of S-field not varying in the x— and y— directions is 
valid as long as spatial variation of S in these two directions is much slower 
than that of perturbation. 



A.l Drift waves associated with density gradients 

Let us impose a small electrostatic potential perturbation at y = as shown 
in Figure A.l with a finite spatial structure in all three spatial directions]^ Fur- 



thermore, consider a case where the dynamics of the perturbation is much slower 
than electron parallel streaming time with its thermal velocity, which is proven 
to be true later. 

In the case that electrons satisfy Boltzmann response, i.e., 5n/n = e^p/Te where 
n {6n) is the mean (perturbed) density, e the proton charge and the electron 



^The lecture was given during the Culham Summer School in 2009. 
^Finite parallel structure is critical for a drift wave to exist 1,26,,164J 



A.l. Drift waves associated with density gradients 




Figure A.l: A schematic of density- 
gradient driven drift waves when 
perturbations in density (black 
solid) and potential (red dashed) are 
in phase. There is no net outward 
particle flux in the x— direction 
while the structures are propagat- 
ing in the +?/— direction (electron 
diamagnetic direction). The mode is 
purely oscillatory, i.e., not unstable. 




Figure A. 2: A schematic of density- 
gradient driven drift waves when 
perturbations in density (black 
solid) and potential (red dashed) 
are not in phase. If the density per- 
turbation leads (lags) the potential, 
then there is net particle flux down 
(up) the density gradient, and the 
growth rate of the mode is positive 
(negative). 



temperature, there is neither net transport of particles in the x— direction nor un- 
stable modes, i.e., the mode is purely oscillatory |164| . This is illustrated with Fig- 
ure A.l The potential structure with a local maximum at y = 0, i.e., ipi > > 0, 
creates a radial electric field which produces E x B drift velocity in the clock- 
wise direction. The particle flux towards positive x— direction (at the top yellow 
box (2)) is r> = [(n> + 5n2) — (n< + ^^2)] ve', while towards negative x— direction 
(at the bottom yellow box 0) is r< = [— (?2> + 5n2) + (n< + ^^2)] ve- Note that ve 
at these two locations is purely in the x— direction. Indeed, we have r> + r< = 0. 

On the other hand, the location of the maximum positive potential is propa- 
gating in the y— direction. Consider a point at y = 5y and constant x denoted as 



in Figure A.l The continuity equation dn/dt + V ■ (nv) = at this location can 



be written as: 

On dn c (f dn 

dt ^ dx 5y B dx^ 
where we have replaced with ve cit y = Sy in the last step. The original per- 
turbed maximum density at y = is Sn = neip/Te, and we can calculate how long 
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A.2. Drift waves associated with ITG and PVG 



it takes to reach this density at the location (2) using Eq. ( A.l ): 

^ c ip dn eip 

6y B dx Te ' 

5y cTe d ^ / * 

Tt = -^9^^"^ = "^- ^^-'^ 

which shows that the potential structure propagates in +?/— direction (electron 
diamagnetic direction) with the speed of Vde- It was mentioned that the drift 
wave dynamics are much slower compared to the electron parallel streaming 
time. This can be validated using Eq. { A.2\ : Vde = ViY^ePe/Ln ^ t'the as pe ^ Ln 
where = | V In ?t,| and pe the electron Larmor radius. 

If the motion of electrons in the parallel direction is impeded by collisions, it 
creates a phase shift between the density and potential structures as shown in Fig- 
ure A.2 In this case, there exists net outward particle transport: the E x B veloci- 



ties at the locations of (T) and (2) are the same as they lie on a constant contour line 
of the potential, but the number of particles at (2) is larger, i.e., more perturbed 
density, than that at ®. Thus, as the structure propagates in the +y— direction, the 
initial perturbation is reinforced, i.e., the mode is unstable. Note that if density 
perturbation lags the potential perturbation, then the mode decays. 



A.2 Drift waves associated with ITG and PVG 

In this section, we describe the ion-temperature-gradient (ITG) and parallel- 
velocity-gradient (PVG) driven modes using the fluid equations of ions where 
electrons are assumed to be in the Boltzmann equilibrium, i.e., we consider ion 
dynamics here. Furthermore, we take the case of a flat density profile (a fi- 



nite density gradient is considered in Appendix A.l) with |VlnTj| = L^^ and 
I V In f/xl = L^^ where Tj is the ion temperature and Uz the incompressible paral- 
lel (f ) plasma flow velocity]^ Again, we consider electrostatic perturbations such 
that 6E = —Vip for which the E x B drift velocity ve = — c (Vip x z) /B. Note 
that i;^ is a perturbed quantity and divergence free by definition. 
The linearized ion continuity equation is 

06n ■ 

— + mViiSUz + Ve ■ Vn, + UzViiSm = 0, (A.3) 
ot 

where we have dropped n{V ve (divergence free), 6niV\\Uz (incompressible) and 
5UzV\\ni (mean density has no gradient in the parallel direction) terms. Here, the 



^Note that Uz is used to denote the 'toroidal' plasma flow velocity in Chapter |5] 
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prefix 6 denotes the perturbed quantity. 

The linearized evolution equation of the parallel momentum|^is 

+ ve ■ VU, + U,V\\6U, j = -\/\\6pi - eniV||(/7, (A.4) 

where we dropped 5UzV\\Uz term. 

The linearized evolution equation of the ion pressure is 

^ + f/.V|)5p. + VE ■ Vp. = -rp,V||(5t/„ (A.5) 

where we have dropped 5UzV\\Pi (no variation of equilibrium pressure in the par- 
allel direction) and —T5piV\\Uz. Here, F is the specific heat ratio. 

The unknown variables in the system are ip, Srii, 5Uz and 5pi. The fourth 
equation is constructed by invoking the Boltzmann electron with quasi-neutrality 
plasma condition: 

5ni ^ Sue ^ n~, (A.6) 

i 

where we assume Ti = T^. This closes the set of equations. 

In our system, a perturbed quantity ^ whose amplitude is taken to be small 
enough can be Fourier decomposed in the y— and z— directions, but such a de- 
composition cannot be done in x— direction as the equilibrium quantities are 
varying in this direction: ^ = ^ (x, t) exp {—iut + ikyU + ikzz) where ^ is the eigen- 
function describing the amplitude of the wave-like structure ^. Thus, using the 
following ansatz, 

^ = (p exp {—iut + ikyy + kzz) , 
Srii = 5ni exp {—iut + ikyU + k^z) , 
5Uz = SUz exp {—iut + ikyy + kzz) , 

5pi = 6pi exp {—iut + ikyy + kzz) , (A.7) 

we look for a condition for a solution of the linearized equations to exist, which 
provides a dispersion relation, with assumptions of ky ^ kz (perturbed quantities 
are elongated along the magnetic field line) and u ^ kzUz, i.e., we make sure that 
kz is small enough such that these assumptions are valid. Note that the E x B 
drift velocity is (we drop the exp term), then, 

Ve = —ikyC^x. (A.8) 
B 



^Perpendicular momentum equation is used to get the E x B drift velocity v e • 
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Rewriting Eqs. (A.3)-(A.6| using Eq. (A.7|, we have 

— u6ni + k^riidUz ~ 0, 
Tiirrii I —UJ6U2 

n 

The dispersion relation we find is 



B 



-Tk^riiTiSUz, 



Ui- 



T- 



(A.9) 



(A.IO) 



where 



ckyTj \ VU^ 
~eB ~C7 



kyPi 



Vthi 



ckyTi 



IVTJ 



eB ( 

i + r- 



i + r)T, 



kyp., 



rrii 



(A.n) 



Using the discriminant A of the cubic function Eq. (A.IO I in to, we can find a 
condition for which an unstable mode exists, i.e., A < such that one root has a 
positive imaginary part: 



A = 4 {klC^^ - k.CsU.uf - 27 [klClu^rY < 0, 

S /I 



:. klCl < k.CsU^u + 



2 \ 2 



(A.12) 



It was mentioned in Appendix |A.l that the density-gradient driven drift wave 
becomes unstable when the density perturbation leads the potential perturba- 
tion, otherwise the wave is purely oscillatory (if they are in phase) or damped 
(if density lags potential). The instability criterion given in Eq. ( A.12[ > can be ex- 
plained in the exact same way. Imaginary part of uj, denoted as tuj, contains the 
phase information for Eq. ( A.9[ >. Thus, if A > 0, we have co[ = meaning that all 
the perturbed quantities are in phase, i.e., purely oscillatory. On the other hand, 
the A < condition (Eq. ( A.12[ >) ensures that there exists a pair of complex conju- 
gate root^so that there exist phase shifts; one of them is growing and the other 
decaying. Note that perturbations in the density and the potential are always in 
^This is alway true for a cubic equation with real coefficients. 
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Figure A.3: A graphical representation of the ITG (x— axis) and PVG (j/— axis) 



instability criterion using Eq. ( A.12) 



phase here via the Boltzmann relation, but the perturbations in the pressure and 
the potential has a finite phase shift if A < 0. 



The dispersion relation Eq. (A.lOl contains three frequencies related to dis- 
tinct physical effects: k^Cg (the sound wave), u^u (convection of sheared flow via 
ve) and cu^t (convection of temperature via ve)- As in any system, a pressure 
perturbation launches a sound wave in our system, but only in the parallel direc- 
tion due to strong magnetic field. Then, we may expect that the system becomes 
unstable if the sound wave cannot keep up the driver of the pressure perturba- 



tion {u^u and/or co^t), and this is exactly what Eq. (A.12 1 dictates. These two 



drivers contribute to destabilize the sound wave together. Figure A.3 illustrates a 
region where unstable modes exist as functions of the sound frequency normal- 
ized driving frequencies. 

We finish this section with a few remarks: (i) both u^u and cu^t require a finite 
ky (see Eq. ( |A.ll l), but Eq. (A.12 1 is not valid for large ky as the finite larmor ra- 
dius effeclj^starts to play a role in which case a proper gyroktnetic equations must 
be used; (ii) shearing rate of the mean perpendicular plasma flow can change the 
instability criterion which was not included in this section; (iii) the smaller the 
kz, the easier to satisfy the instability criterion, and the smallest possible /c^ in a 
tokamak is the inverse of the connection length. 



^If scales of perturbed potential is smaller than the ion Larmor radius, then the gyro-motion 
effectively averages out the perturbed potential. 
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Examples of a correlation function and a 

power spectrum 



In this chapter, we describe, step by step, how we generate a correlation func- 
tion and a power spectrum in this work from the raw BES data. 

B.l Generating a correlation function 

(1) Get the raw BES data and frequency-filter it from 20 to 100 kHz as shown in 



Figure B.l 'a) and (b). 

(2) Select the time range where correlation functions are to be generated. In Fig- 
ure 



B.l Id), we selected time range of 0.180 — 0.185 s. 
(3) Divide the selected time range into subwindows. In this example, we divided 
the time range into 50 subwindows with 100 /is duration for each sub window. 



(4) Using Eq. (5.9), generate a correlation function for each subwindow. In our 



example, we, then, have 50 correlation functions as shown in Figure B.2[ a) with 



black lines. Note that we have cross-correlation functions using two poloidally 
separated channels, i.e., Ch. 10 and Ch. 18 whose separation distance is 2.0 cm. 
(5) Average the 50 correlation functions to obtain the averaged correlation func- 



tion shown as red circles in Figure B.2 'a). 

(6) Uncertainty of the averaged correlation function is estimated as the standard 
deviation divided by the square root of number of subwindows II165L i.e., VSO in 
our example. 



B.2. Generating a power spectrum 
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Figure B.l: (a)-l and (a)-2 show raw BES data of Ch. 10 and Ch. 18 from shot 
#27267; (b)-l and (b)-2 show frequency filtered ([20, 100] kHz) BES data of (a) 
from 0.180 to 0.185 s, i.e., red vertical dashes in (a). Note that Ch. 10 and Ch. 18 
has a poloidal separation distance of 2.0 cm with no radial separation. 



B.2 Generating a power spectrum 



(1) Get the raw BES data as shown in Figure B.l 'a) 



(2) Select the time range where power spectra are to be generated. In Figure 



B.l 'b), we selected time range of 0.180 — 0.185 s. 

(3) Divide the selected time range into subwindows. In this example, we divided 
the time range into 50 subwindows with 100 jjs duration for each subwindow. 

(4) Generate a cross-spectrum for each subwindow as = 
FT {fi{t)} FT* {gi{t)}, where FT {} and FT* {} are the Fourier-transform 
operator and its complex conjugate, respectively. Here, the subscript i denotes 
the ith subwindow, and / (t) and g (t) for raw BES data from 0. 180 to 0. 185 s in our 



example. Black lines in Figure B.2p 3) shows the magnitude of 50 cross-spectra. 



i.e., 1 5'/^ I, where \S\ returns the magnitude of the S. 

(5) Average the 50 15/^1 to get the averaged power spectrum shown as red 
circles in Figure B.2 'b). (Differences between red circles and green triangles are 
explained later.) 

(6) Uncertainty of the averaged power spectrum is estimated as the standard 



deviation divided by the square root of number of subwindows 11651 - 
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Figure B.2: (a) Cross-correlation functions and (b) cross-power spectra between 
Ch. 10 and Ch. 18 shown in Figure B.l 'b). 



We note that red circles in Figure B.2 h) are calculated as the average of the \Sl'\ 
while green triangles are calculated as the magnitude of the averaged S{^. These 
two estimates give different results because taking the magnitude of complex 
number is not a linear operator. Following Bendat and Piersol |165| , we actually 
take the magnitude of the averaged in this work, i.e., green triangles, as the 
power spectrum of the signal. The reason we show 15/^1 (black lines) and the 
averaged | S^^ \ (red circles) is because of the difficulty of plotting for it being 
complex numbers. Note that the uncertainty for the magnitude of the averaged 
gf9 (green triangles) is estimated by first dividing the standard deviations of real 
and imaginary parts of the by the square root of the number of subwindows, 
then by propagating these uncertainties to get the corresponding uncertainty 
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Definitions 



Fields 



B Magnetic field 

Toroidal component of the magnetic field 

Bg, Bp Poloidal component of the magnetic field 

E Electric field 

J Current density 

(f Perturbed electrostatic potential 

Lengths 

pi Ion Larmor radius = fthi/^i 

Pe Electron Larmor radius = Vthe/^e 

P* Pi/ a 

Ln Density gradient scale length = | V lnn|~^ 

Lt Temperature gradient scale length = | V In T|~^ 

Lt. Ion temperature gradient scale length 

Lt^ Electron temperature gradient scale length 

L* Minimum of L^, and Ln 

Lu Flow velocity gradient scale length = | V In f/|~^ 

\x Radial correlation length of Gaussian "eddies" in Chapters |4] and |5] 

Ay Poloidal correlation length of Gaussian "eddies" in Chapters |4] and |5] 

fcx Perpendicular wavenumber of plasma turbulence 

fc| Parallel wavenumber of plasma turbulence 

£x Perpendicular correlation length of plasma turbulence 

£|l Parallel correlation length of plasma turbulence 

^xr (^-R Radial correlation length of plasma turbulence 

£y Perpendicular (on a flux surface) correlation length of plasma turbulence 



Times (frequencies) 



Iz Poloidal (vertical) correlation length of plasma turbulence 

/c^ I 

ky '2i7T ^ £y 

A The connection length = nrB/Bp 



Times (frequencies) 



te Energy confinement time 

Hi Ion gyro-frequency = eB/niiC 

f2e Electron gyro-frequency = eB/nieC 

Tiife Lifetime of Gaussian "eddies" in the moving frame in Chapters |4] and |5] 

^poak ■ ■ ■ Time delay at which the cross-correlation function reaches its maximum 

Tc Correlation time of plasma turbulence 

Nonlinear time of plasma turbulence 

r^^ Nonlinear time associated with the non-zonal component of plasma 

turbulence 

Tst Particle (ion) streaming time along the parallel direction = A/fthi 

Tsh Perpendicular velocity shear time = [{Bp/ B) dU^/dr]~^ 

tm Magnetic drift time (ix/Pi) R/vtu 

T^i Drift time associated with ion- temperature gradient (iy/pi) LTjvtu 

r*e Drift time associated with electron-temperature gradient {ly/ pe) Lt^ /vthe 

T^n Drift time associated with density gradient = {ty/ pi) Ln/v^-<a.i 

r* Minimum of r*j and r*„ 

Uii Ion collision frequency 

u^i Normalised ion collision frequency = uuTst 

uj^u Frequency of the drift wave associated with the parallel flow shear 

ijj^T Frequency of the drift wave associated with the temperature gradient 

7£; Perpendicular flow shear = {Bq/B) dU^/dr 

Normalised perpendicular flow shear = 'jETst 



Velocities 

c The speed of light 

Vthi Ion thermal velocity = ^J^Tijmi 

f the Electron thermal velocity = ^J'lT^|m^ 
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Other 



Cs Sound speed = (1 + F) Ti/rrii 

U Plasma flow velocity 

f/|l Parallel component of plasma flow velocity 

U± Perpendicular component of plasma flow velocity 

f/^ Toroidal component of plasma flow velocity 

Uz Toroidal component of plasma flow velocity in Chapter |5] 

Uy Poloidal component of plasma flow velocity in Chapter |5] 

ve, Su± Fluctuating E x B drift velocity 

vd Magnetic drift velocity 

Vy^^ Apparent poloidal velocity of measured fluctuating density patterns 

Vy Apparent poloidal velocity of Gaussian "eddies" in Chapters 

Vy^^^ RMS value of fluctuating Vy of Gaussian "eddies" in Chapters 



4 


and 


5 


4 


and 


5 



Other 



n Plasma density 

6n Fluctuating plasma density 

rij Ion density 

Ue Electron density 

T Temperature ineV 

Ti Ion temperature ineV 

Te Electron temperature in eV^ 

p Plasma pressure 

rui Ion mass 

mg Electron mass 

/3 The ratio of the plasma pressure to the magnetic field energy density 

e Proton charge 

/ Mean part of the detected photon intensity by the 2D BES 

61 Fluctuating part of the detected photon intensity by the 2D BES 

Imhd RMS value of MHD modes in synthetic data in Chapters [1] and [s] 

r Specific heat ratio 

Qturb Gyro-Bohm-normalised turbulent ion heat flux 
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Geometrical quantities 



Geometrical quantities 



ijj Flux surface label 

R Distance from the centre of a torus (major radius) 

r Distance from the magnetic axis of a flux surface (minor radius) 

Z Height from the midplane of flux surfaces 

Half diameter of a flux surface ip at the midplane 

a of the LCFS, a measure of total plasma size 

i?o ^ at the point where — s> 

Eq Inverse aspect ratio of a tokamak, i.e., = a/Ro 

e Local inverse aspect ratio, i.e., e = r/R 

q The safe factor ~ §^ for large aspect ratio tokamak 

a Pitch angle of local magnetic field, tan a = 

Acronyms 

BBS Beam Emission Spectroscopy 

CCFE Culham Centre for Fusion Energy 

CCTD Cross-correlation Time Delay 

CUDA Compute Unified Device Architecture 

CXRS Charge eXchange Recombination Spectroscopy 

ETG Electron Temperature Gradient 

GPU Graphical Processing Unit 

ITB Internal Transport Barrier 

ITG Ion Temperature Gradient 

LCFS Last Closed Flux Surface 

LoS Line-of -Sight 

MAST Mega Amp Spherical Tokamak 

MHD Magnetohydrodynamic 

MSE Motional Stark Effect 

NBI Neutral Beam Injection 

PSF Point-Spread-Function 

PVG Parallel Velocity Gradient 

RMS Root-Mean-Square 

START Small Tight Aspect Ration Tokamak 

TEM Trapped Electron Mode 
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